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' baseline  rocket  and  parametric  variations  thereof.  These  results  indicate 
that  transverse  flexibility  contributes  significantly  to  rocket  mallaunch 
when  the  spin  rate  is  high,  since  the  long-period  part  of  the  angular  rate 
at  launch  is  directly  proportional  to  spin  rate  and  to  the  magnitude  of 
transverse  deformation.  Effects  of  thrust  misalignment  are  considered  in 
comparing  results  obtained  from  solutions  to  the  linear  equations  with  flight 
data  for  a particular  rocket.  Equations  for  motion  of  the  two-body  rocket  on 
a "rigid"  (non-moving)  launcher  are  also  derived.  Results  obtained  by  numeri- 
cally integrating  these  equations,  are  presented.  The  second  and  more  complex 
model  is  that  of  a continuous,  non-uniform,  slender  rod  subject  to  "internal" 
thrust.  Equations  for  motion  on  and  off  a "rigid"  launcher  are  derived  and 
converted  to  ordinary  differential  equations  by  vising  plnned-plnned-free 
mode  shapes  for  the  on-*launch  phase  and  free-free  mode  shapes  for  the  free- 
flight  phase.  Conversion  from  the  on- launcher  equations  to  the  free-flight 
equations  is  explained  and  preliminary  results  are  discussed. 
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PREFACE 


This  report  contains  a description  of  the  work  accomplished  under  U.S. 
Army  Contract  DAAK40-C-77-0125,  during  the  period  13  June  1977  to  30  Sep- 
tember 1978,  for  the  U.S.  Army  Missile  Research  and  Development  Command, 
Redstone  Arsenal,  Alabama.  Mr.  Dean  E.  Christensen,  the  Contracting 
Officer's  Technical  Representative,  was  of  great  assistance  during  this 
effort,  through  both  his  encouragement  and  technical  suggestions. 

The  authors,  other  than  the  first,  contributed  to  the  contractual 
effort  in  the  following  ways.  Mr.  Carl  Evans  wrote  a major  segment  of  the 
computer  code  for  the  continuous  rocket  model  simulation  and  also  worked 
out  some  of  the  mathematical  details  of  the  model  before  he  had  to  enter 
the  U.S.  Air  Force  in  the  spring  of  1978.  Mr.  James  B.  McCurry  assisted 
in  programming  codes  for  the  linearized  two-body  rocket  model,  as  well  as 
providing  help  on  working  out  the  two-body  rocket  model  equations  for  motion 
on  a rigid  launcher.  Others  who  made  contributions  are  Mr.  Steve  Lowe, 
undergraduate  research  assistant,  who  performed  exceptionally  well  as  a 
programmer,  Mr.  Grant  Castleberry  and  Mr.  James  Thompson,  graduate  research 
assistants,  who  also  assisted  in  programming,  and  Mr.  Gene  Holloway,  graduate 
research  assistant  who  drew  most  of  the  figures.  Finally,  Mrs.  Marjorie 
McGee  deserves  special  recognition  for  her  excellent  typing  of  the  manuscript. 


John  E.  Cochran,  Jr. 
Project  Leader 
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SECTION  1.  INTRODUCTION 


Recently  studies  have  been  initiated1  with  the  goal  of  improving  the 
accuracy  of  free-f light  rockets.  Although  significant  reductions  in  disper- 
sion of  free-flight  rockets  have  been  accomplished  by  utilizing  modern 
manufacturing  methods,  by  using  various  techniques,  such  as  spinning  the 
rocket  to  minimize  the  effects  of  mechanical  thrust  misalignment , and  by  using 
launchers  which  provide  release  of  the  rocket  without  producing  appreciable 
tipoff,  further  improvement  in  free-flight  rocket  accuracy  will  forseeably 
require  more  detailed  analysis  of  factors  which  were  once  considered  relatively 
insignificant.  One  such  factor  is  that  of  transverse  flexibility  of  the  rocket. 

Since  no  structure  is  perfectly  rigid,  all  rockets  are  flexible  to  some 
extent.  In  the  past,  however,  most  free-flight  rockets  were  designed  with 
length  to  diameter  (L/D)  ratios  less  than,  say,  fifteen.  Furthermore,  the 
structures  of  these  rockets  tended  to  be  dense.  Such  rockets  are  so 
relatively  rigid  that  they  may  be  modeled  well  as  variable-mass  rigid  bodies.2 
Advances  in  materials  science  have,  however,  presented  the  free-rocket 
designer  with  strong, light-weight  alloys  and  composite  materials.  Free-flight 
rockets  constructed  using  such  materials  will,  in  general,  be  lighter 
(warhead  excluded)  than  their  predecessors  and  often  more  flexible.  Further- 
more, efficiency  in  terms  of  aerodynamic  design  dictates  large  L/D  rockets 
which  are  now  structurally  possible  configurations,  but  which  tend  to  be 
more  flexible  than  the  previous  rockets  with  smaller  L/D's. 

The  magnitude  of  the  effect  of  the  flexibility  of  a spinning  free  rocket 
on  dispersion  has  not  been  established  accurately.  In  fact,  although  a good 
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deal  of  work  has  been  expended  in  the  area  of  determining  the  effects  of 
flexibility  on  the  stability  and  control  of  rockets  of  the  launch  vehicle 
type,3,4’5  relatively  little  work  concerning  flexible,  spinning  rockets 
has  been  done.  In  this  regard,  an  article  by  Reis  and  Sundberg6  addresses 
the  problem  of  aeroelastic  bending  of  spinning  sounding  rockets  using  a simple 
two-rigid-body  model  similar  to  one  of  these  discussed  in  the  body  of  this 
report,  while  Womack,  et.  al . 7 discuss  the  use  of  a more  complex  model  based 
on  the  use  of  normal  modes  for  modeling  the  flexibility  of  such  rockets.  A 
general  approach  to  the  flexible  spinning  rocket  problem  is  taken  by  Meriovitch.8 

The  problem  treated  here  is  different  from  the  problem  of  modeling  the 
dynamics  of  a flexible,  spinning,  sounding  rocket,  because  in  the  launch  of  a 
free-f light  rocket,  the  rocket  passes  from  a constrained  condition  to  a free 
condition.  This  transition  is  of  utmost  importance  for  free  rockets,  because 
their  accuracy  is  greatly  affected  by  their  state  at  the  time  of  release 
from  the  rocket  launcher. 

The  primary  purpose  of  the  effort  reported  on  here  is  to  provide  a basis 
for  determining  the  effects  of  transverse  flexibility  of  free-flight  rockets 
on  their  motion  during  the  launch  phase  and  during  the  free-flight  phase  up  to 
the  time  when  aerodynamic  reactions  become  significant.  A two-pronged 
approach  was  adopted.  On  one  hand,  the  most  simple  model  which  it  was  felt 
would  give  meaningful  results  was  formulated  and  results  obtained  therefrom. 
Concurrently,  a fairly  sophisticated  model,  based  partly  on  the  results  of 
Ref.  9 was  developed. 

The  simple  model  of  a flexible  rocket  consists  of  two  rigid  bodies 
coupled  in  such  a way  that  only  transverse  relative  rotation  is  allowed  and 
that  rotation  is  resisted  by  linear  torsional  springs  (with  viscous  damping 
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also  present).  This  model  Is  described  in  Section  2.  In  that  section, 
nonlinear  equations  of  motion  for  the  free-flight  phase  are  derived  and 
linear,  constant-coefficient  approximations  to  these  obtained.  Results 
obtained  by  solving  the  linear  equations  in  closed  form  are  presented 
also.  Furthermore,  nonlinear  equations  of  motion  for  launch  from  a non- 
moving, or  "rigid,"  launcher  are  derived  and  some  results  of  numerical 
solutions  to  these  are  presented  and  discussed. 

The  more  complex  model  is  described  in  Section  3.  This  model  is  that 
of  a slender,  elastic  body  which  has  acting  on  it  a thrust  force  and 
gravity  (No  aerodynamic  reactions  are  considered.).  Hybrid  differential 
equations,  i.e.,  differential  equations  which  contain  both  partial  and 
ordinary  derivatives,  are  derived.  These  are  replaced,  approximately,  with 
ordinary  differential  equations  by  using  the  method  of  assumed  modes.10 
Two  different  types  of  mode  shapes  are  used.  For  the  launch  phase,  pinned- 
pinned-free  mode  shapes  for  a uniform,  slender  beam  are  used  because  of 
the  desire  to  model  the  launcher  (considered  "rigid")  constraints.  For  the 
free-flight  phase,  free-free  mode  shapes,  also  for  a uniform,  slender  beam, 
are  used.  Furthermore,  since  instantaneous  principal  axes8  are  used  for  the 
free-flight  phase  as  well  as  the  different  mode  shapes,  the  equations  needed 
to  transform  the  "final  conditions"  of  the  launch  phase  into  "initial  con- 
ditions" for  the  free-flight  phase  are  given  also.  Due  to  the  fact  that  the 
digital  computer  code  which  is  used  to  produce  numerical  solutions  in  both 
sets  of  equations  was  not  completely  debugged  by  the  due  date  of  this  report 
no  results  obtained  from  it  are  presented.  However,  some  comments  are  made 
in  regard  to  "anticipated"  results. 

In  Section  4,  conclusions  reached  during  the  course  of  this  investiga- 
tion are  stated  and  recommendations  for  further  study  are  made. 


SECTION  2.  TWO-BODY  FLEXIBLE  ROCKET  MODEL 


2.1  Description  of  the  Model 

The  two-body  physical  model  of  a free-f light  rocket  is  shown  in  Fig.  1. 

The  model  is  composed  of  two  almost  axisymmetric  rigid  bodies.  These  bodies 
are  connected  at  point  R in  such  a manner  that  the  forward  body  (hereinafter 
called  Body  2)  may  rotate  relative  to  the  aft  body  (hereinafter  called  Body  1) 
transversely  with  respect  to  the  symmetry  (if  no  dynamic  imbalance  is  present) 
axis  of  Body  1.  This  relative  rotation  is  elastically  restrained  and  viscous 
damping  is  also  present. 

The  point  P is  located  at  the  aft  end  of  Body  1 and  the  center  of  mass  of 
Body  1,  denoted  C^>  is  located  with  respect  to  P by  the  vector  r^.  The  vector 
connect;,  the  points  P and  R and  the  vectors  r^  and  are  directed  from  R to 
C2  (the  center  of  mass  of  Body  2)  and  from  to  Q,  respectively.  The  points  P 
and  Q correspond  to  the  points  at  which  the  rocket  is  supported  while  it  is  on 
the  launcher. 

During  the  launch  phase,  it  is  assumed  that  the  launcher  does  not  move. 
Hence,  it  is  referred  to  as  a "rigid"  launcher.  Thus,  while  on  the  launcher, 
the  rocket  is  constrained  so  that  the  points  P and  Q translate  along  a fixed 
straight  line.  In  addition  to  the  constraint  forces  acting  at  P and  Q,  the 
rocket  is  subjected  to  external  forces  due  to  thrust  and  gravity  while  it  is  on 
the  launcher.  However,  its  mass  is  assumed  constant  during  this  phase. 

During  the  free-f light  phase,  a free-f light  rocket  is  generally  acted  on 
by  gravity,  aerodynamic  reactions  and  forces  and  moments  due  to  the  flow  internal 
to  the  rocket.  At  present,  it  is  desired  to  maintain  simplicity  in  our  model. 
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Hence,  the  aerodynamic  reactions  as  well  as  all  the  forces  and  moments  on  the 
rocket  due  to  internal  flow,  except  the  thrust,  are  neglected.  These  should 
be  relatively  small  compared  to  the  thrust  and  weight  during  the  first  two  or 
three  tenths  of  a second  of  flight.  The  thrust  force  is  assumed  to  act  only 
on  Body  1 and  its  magnitude  and  direction  relative  to  the  longitudinal  axis  of 
Body  1 are  assumed  constant.  Furthermore,  the  mass  of  the  rocket  is  assumed 
constant  during  this  phase  as  in  the  launch  phase.  The  assumption  of  constant 
mass  should  be  fairly  valid  for  the  time  period  immediately  after  launch  during 
which  optical  lever  data  is  acquired. 


Figure  1.  Simple  two-body  rocket  model. 
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2.2  Nonlinear  Equations  for  the  Free-Flight  Phase 

The  equations  of  motion  for  the  physical  model  described  for  the  off-the- 
launcher,  or  free-flight  phase, + may  be  obtained  by  applying  Newton's  laws  of 
motion  directly,  or  by  using  an  energy-based  method  such  as  Lagrange's.  Here 
we  choose  the  direct  approach  in  the  sense  of  Ref.  2.  Equations  of  motion 
for  the  system  of  two  bodies  are  first  obtained.  Next,  the  equations  of  motion 
for  Body  2 are  derived. 

Explanation  of  Notation 

A dextral,  orthogonal,  coordinate  system,  OXYZ,  fixed  to  the  earth's 
surface  is  used  as  the  reference  for  all  motion.  A Body  1-fixed  coordinate 
system,  311,1  3 corresponding  Body  2-fixed  system,  c2X2y2Z2’  are  defined 

such  that  the  x^-axis  is  the  longitudinal  axis  of  Body  j.  The  y^-  and  z^-axes 
are  arbitrarily  oriented  and  the  y^~  and  z2~axes  lie  in  the  x^y^-  and  x^z^-planes, 
respectively,  when  the  x^-  and  x2~axes  are  colllnear.  The  system  Cxyz  is  defined 
such  that,  although  its  origin  is  the  center  of  mass  of  Bodies  1 and  2,  its  axes 
are  aligned  with  the  C^y^  system  (see  Fig.  2).  The  vectors  ^ and  £2  shown 
in  Fig.  2 are  defined  by  the  two  equations. 


£2  ~ -2  * -1  + -2 

(1) 

- 

and 

1 

ml  £l  + m2  -2  “ - * 

(2) 

where  m^  is  the  mass  of  Body  j.  Thus, 

£l  " + £2^ 

(3a) 

+The  free-flight  phase  is  considered 

first  in  this  report  because  it  was  ^ 

considered  first  in  time  during  this  contractual  effort. 
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In  addition  to  the  notation  thus  far  defined,  we  shall  use  a dot  over  a 
matrix,  or  a vector,  to  denote  total  differentiation  of  that  quantity  with 
respect  to  time,  while  a small  circle  over  a vector  will  be  used  to  denote 
differentiation  - with  respect  to  time  of  only  the  components  of 
that  vector  as  resolved  in  the  Cxyz  system.  Furthermore,  the  notation  6V/6t 
is  used  to  denote  the  derivative  of  the  components  of  a vector  quantity  V as 
resolved  in  the  C2X£y2Z2  coor<**-nate  frame*  A bar  (line)  under  a quantity  is 
used  to  indicate  that  such  quantity  is  either  a vector  or  a column  matrix,  while 
two  bars  (lines)  under  a quantity  indicate  that  it  is  a dyadic  or  a square  matrix. 
Finally,  a tilde  above  a 3x1  column  matrix  is  used  to  indicate  that  a skew- 
symmetric  3x3  matrix  is  to  be  formed  from  its  elements  and  zeroes.  For  example, 
given 

(4) 

(5) 

Equations  of  Motion  for  the  System 

Newton's  laws  of  motion  for  a particle  may  be  used  to  show  that  the  time  rate 
of  change  of  the  angular  momentum  of  the  system  composed  of  Bodies  1 and  2 about 
C is  equal  to  Tc>  the  external  torque  about  C.  The  angular  momentum  of  the  system 
about  C,  H , may  be  expressed  in  the  form, 

l 


0 A 


0 

o. 


-n. 


* 


Sc  = il  * - + ^2  *(-  + -2/l)  +mi  P-1*P-1 


+ m2  P2  * P2  • 

By  using  Eqs.  (3),  we  may  rewrite  H as 

-c 

Sc  ■ £1  * 2 + *2  *(-  + -2!\}  + a(-l  + H2)x(^1  + (7) 

where  a ■ m^m^/M.  Since  Hc*Tc>  it  follows  from  Eq.  (7)  that 

T - I , • Q + Q x I , • a + I.  ‘(fiQ/fit  + fiu  . /fit) 

+ (n  + “2/1^  x ^2  + -2/1^  + °^-l  + -2*  X ^-1  + ^2^  ’ <8) 

The  right-hand  side  of  Eq.  (8)  is  actually  more  complicated  than  it  might  first 
appear,  because  and  r2  are  rather  complicated  expressions  when  expanded  out 
in  full. 

It  is  convenient  to  convert  Eq.  (8)  into  its  matrix  counterpart  at  this 
point.  If  we  use  the  unit  vectors  i»  5 and  £ (which  are  parallel  to  the  x- , 
y-  and  z-axes,  respectively)  as  our  basis  vectors  (Body  1 basis)  for  writing  the 
matrix  equation  and  want  to  express  1 2 and  r^  in  the  Body  2 basis  (unit 

A A A 

vectors  12»  J2  and  k2),we  need  a transformation  matrix  to  relate  the  two  bases. 

A suitable  matrix  is 


C0,C0, 

C0.S0 

2 3 

2 

-80, 

C0  , 

3 

3 

C0„S0„ 

-S0,S0 

3 2 

2 

N 


where  6^  and  6^  (see  Fig.  3)  are  angles  of  relative  rotation  of  Body  2 with 


respect  to  Body  1 about  the  y - and  z -axes,  respectively,  and  c0-  cos  0 and 


s0«  sin  0.  The  matrix  A thus  relates  the  unit  vector  triads  (i,1,k)  and 


(^2 ’-2-2^  throu8h  the  equation. 


Figure  3.  Relative  orientation  of  bodies  1 and  2 
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It  can  be  shown  that 


A * “ —2/1  * * 


and,  since  (a  number  under  a matrix  denotes  the  basis  used) 

£ - A £ , 

2 1 

the  matrix  counterpart  of  the  vector  <5ft/<5t  is  ”“2/1  ■ ^ + & 5 • This  fact  may 
be  used  to  get  the  results,* 

t - i ft  + 5 1 ft  + at  i,  (a  ft  a>, a c + <L0/.) 


-c  * - 


1,  Wj  i, 
2 " 2 ' 


+ (^  AT  + AT  “2/1^  ^2  n + -2/1^ 

+ o^  + AT  r2  A)  t-  2 + 2 2 I (13) 

- iT  §2  2 - “2/1  i 2 + “2/l> 

+ (ft  + A^  §2/1  A)  (2  + A^  <^2/1  «) A^ 
where  the  Body  1 basis  is  used  unless  otherwise  indicated.  Also 

T - - [p(L  + AT  r.  A)  + r.]F  (14) 

I " 2 I 

where  F^,  is  the  thrust. 

Equation  (13)  governs  the  motion  of  the  system  of  two  bodies  about  the  system 
center  of  mass.  It  can  be  considered  an  equation  from  which  ft  may  be  obtained. 

I 

An  additional,  independent, matrix,  differential  equation  for  ^2/1  is  needed. 

*A  superscript  T is  used  to  denote  the  transpose  of  a matrix. 
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Relative  Motion  of  Body  2 


A differential  equation  governing  the  motion  of  Body  2 with  respect  to 

Body  1 may  be  obtained  by  using  the  fact  that  the  time  rate  of  change  of  the 

angular  momentum  of  Body  2 about  R,  minus  the  mass  of  Body  2 times  the  cross 

product  of  r^  and  the  acceleration  of  the  point  R,is  equal  to  the  torque  about 

R.  To  see  that  this  is  the  case,  we  write  a,  , the  acceleration  of  dm- 

-anu  z 

(an  element  of  mass  of  Body  2)  in  the  form. 


2d " 5P  + £i  + *i  + l2  + £2  * (] 

where  p£  is  a generic  position  vector  from  to  dm^.  By  forming  the  cross 
product  of  r ^ + P*£  with  each  side  of  Eq.  (15)  and  integrating  over  the  mass 
of  Body  2,  we  get+ 


m2  E2  x + Ii  + ij)  + m2  -2  X E2 


+ J * Pj  dm2  “ / (£2  + £2?  X -dm  dm2 

m2  m2  2 


/ (rj  + £^)  x a^  dm2  ■ T , the  torque  on  Body  2 about  R,  (17) 

m2  m2 

and  the  sum  of  the  second  and  third  terms  on  the  left-hand  side  of  Eq.  (16)  is 
the  time  rate  of  change  of  the  angular  momentum  of  Body  2 due  to  its  rotation 
about  the  point  R.  Also,  the  acceleration  of  R is  R + r,  + i. . 

-p  -1  -1 

Now,  we  first  introduce  the  centroidal  inertia  dyadic  of  Body  2 and  the 
equation  for  translational  motion  of  m^t 

+Note  that  / p*  dm2  2 2»  virtue  of  the  definition  of  p*. 
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We  can  rewrite  Eq.  (20)  in  the  form, 


-21  - + -22  -2/1  “ -2  * 


(21) 


where  the  definitions  of  anc*  T2  ^°^^-ow  from  Eq.  (19).  Then* 


-2/1  " -22  -2  ' -22  -22  2 ‘ 


In  a similar  manner,  Eq.  (13)  can  be  rewritten  in  the  form. 


-11  - + -12  -2/1  " -1  * 

Thus,  by  using  Eq.  (22)  to  eliminate  <£2/1  from  Eq‘  (2^)  » we  8et 
[-ll  ~ -12  -22  -21]  2 " -1  " -12  -22  -2  ’ 


(22) 


(23) 


(24) 


Similarly,  by  solving  Eq.  (23)  for  0 and  using  that  solution  in  Eq.  (21),  we 
obtain 


[-22  " -22  -11  -12 ] -2/1  " -2  ” -22  -11  ^1  * 


(25) 


The  derivatives  0 and  £2/1  can  be  obtained  from  Eqs.  (24)  and  (25)  by  simply 
inverting  the  3x3  matrices  in  square  brackets  and  premultiplying  each  equation 
by  the  proper  invei  se  matrix. 

Kinematic  Equations  for  Rotation 

The  orientation  of  the  Cxyz  coordinate  frame  is  defined  by  the  Euler  angles 
'f , 0,  and  $.  As  shown  in  Fiy.  4,  a 3-2-1  rotation  sequence  is  used.  The 
kinematic  equations  for  these  angles  are 


'f  - (flj  s*n  * + by  cos  *)/cos 


(26a) 


A superscript  -1  is  used  to  denote  the  inverse  of  a square  matrix. 


0 - n2  cos  4>  - ft3  Sin  4> 

and 

1 - + (n2  Sin  $ + «3  COS  *)tan  0, 

where  the  Q.t  j-1,2,3,  are  the  elements  of  Q. 

1 1 


(26c) 


Figure  4.  Orientation  of  body  1. 


The  orientation  of  Body  2 relative  to  Body  1 is  defined  by  the  previously 
introduced  angles  02  and  63.  These  angles  are  used  in  a 3-2  rotation  sequence, 
as  shown  in  Fig.  3.  From  the  geometry  on  Fig.  3,  one  may  deduce  the  kinematic 


equations. 
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03  - w^/ cos  02 


(27a) 


2 * 


(27b) 


from  which  these  angles  may  be  determined.  Here,  and  are  elements  of 

T 

u)2y^>  i.e. , 3 (^  “2  • 

Translational  Motion 

Although  it  is  of  secondary  importance  in  this  investigation,  for  the  sake 
of  completeness,  we  shall  set  out  differential  equations  for  the  translational 
motion  of  the  two-body  rocket's  center  of  mass.  In  point  of  fact,  the  equation 
governing  translational  motions  of  C is  Eq.  (18).  However,  we  chose  to  write 
it  in  the  less  complicated  vector  form, 

V =■  - OxV+F/M  , (28) 

where  V is  the  velocity  of  C.  A matrix  form  of  Eq.  (28) (Body  1 basis)  is 


v * a v 

+ eg  + ft/m 

9 

(29) 

where  g * 

(0 

T 

0 g)\ 

Also,  in  Eq.  (29) 

the  matrix  C 

s 

is  a transformation 

matrix  which 

relates  the 

fixed  OXYZ  basis 

and  the 

Body  1 

basis.  Explicitly, 

r 

cQcV 

cQa'V 

— S0  “ 

C * 

-C't’S'i'+S'l’C'i'S© 

ctcH'+s'ts't'sG 

S$C0 

• 

(30) 

_s$s1'+c<l>c'i'  sin9 

-s<f>c'H-ei>s4's0 

C$C0. 

Furthermore,  V = (u  v w)  . Notice  that,  although  the  values  of  R^,  r^, 
• • 

and  r2  at  the  instant  of  release  of  the  rocket  from  the  launcher  are 
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(32) 


for  the  transformation  matrix  A and  also  may  approximate  Fq.  (27a)  with 

®3  2 “3  • (33) 

Our  second  assumption  Is  that  the  angles  0 and  4*  are  also  small  enough  that 
the  sines  of  these  angles  are  approximately  equal  to  the  ingles  in  radians  and 
that  the  cosines  of  the  angles  are  each  approximately  unity.  The  assumption 
that  9 is  small  limits  the  launcher  elevation  to  small  angles.  However,  at  small 
values  of  0,  the  effects  of  bending  of  the  rocket  due  to  its  weight  will  be 
greatest,  so  that  this  is  a conservative  assumption. 

With  small  0 and  f,  we  get  the  following  approximation  to  C: 


“ 

i 

Y 

-0‘ 

c * 

-¥c* 

+ 

0s4> 

s$ 

• 

(34) 

■e- 

0) 

l 

+ 

0C$ 

-8$ 

c<t>. 

1 

We  also  may  simplify  Eqs.  (26)  to  get 


v * n2  s«  + n3  c<t. 


(35a) 


0 =■  Jlj  c$  - s$ 

and 

1 * + i 0 . 


(35b) 

(35c) 


Now,  if  V and  0 are  to  remain  small,  f and  0 must  be  small.  Furthermore, 
• • 

for  a "good"  launch,  V and  0 are  small  compared  to  one  radian  per  second.  It 


r 


i 
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then  follows,  from  Eqs.  (35),  that  ^ and  0-j  are  the  same  magnitude  as 
• • 

6 and  V and  hence  are  "small." 

We  shall  further  assume  that  >>  4*  0.  Thus,  Eq.  (35c)  can  be  further 
simplified  to  get  the  approximate  equation, 

♦ ■ nx  . (36) 

• • 

The  assumption  that  4*  and  0 are  small  compared  to  unity  makes  it  possible 
for  us  to  simplify  the  equations  of  rotational  motion  significantly.  First  of 
all,  we  can  discard  any  products  of  0^  an<*  ^ ’ anc*  also  any  squares  of  these 
variables,  from  the  equations.  Furthermore,  it  generally  is  possible  to  assume 
that  and  are  also  less  than  unity,  although  for  large-amplitude,  high- 
frequency  vibration  they  may  have  magnitudes  near  unity  at  times.  *"  Thus,  any 
products  of  t>»2  an<^  ^3  or  of  b>2  and  and  or  may  be  discarded  as  being 
higher  order. 

The  terms  which  contain  0^  and  0^  are  retained  for  the  present.  These  may 
be  very  large  compared  to  unity  if  the  rocket  is  spinning. 

Although  more  generality  could  be  maintained  with  additional  effort,  it  is 
assumed  here  that  the  following  geometric  equations  hold: 


Hi  ' (ri 

0 

or  , 

(37a) 

Hi  ’ <li 

0 

0)T 

(37b) 

§2  ' <r2 

0 

0)T  , 

(37c) 

and 


where  r^,  and  r^  are  constants.  These  equations  imply  that  the  centers  of 
mass  of  the  bodies  both  lie  on  the  x-axis  when  the  rocket  is  not  "bent."  ■ 


Generally,  lower  amplitudes  are  associated  with  higher  frequencies. 


r 


■yr,' 


The  inertia  matrices  are  assumed  to  have  the  forms 


which  of  course  imply  that  the  bodies  are  axisymmetric  and  dynamically  balanced 
The  torque  T2/1  is  considered  to  be  a linear  function  of  «2>  0^  02  and  6 and 

|t  " ^FT  azFT  ~ayFT^  ’ where  ay  and  az  are  mechanical  thrust  misalignmen 

angles  as  shown  in  Fig.  5. 


Figure  5.  Mechanical  thrust  misalignment  angles 


After  a good  deal  of  algebra  and  the  use  of  the  above  stated  assumptions,  the 
equations  of  rotational  motion  may  be  reduced  to  the  forms. 
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I 


I x « B x + T, 


where  T is  the  x-component  of  the  torque  T and  I “A  + B,  is  the  moment  of 
s -c  A 1 1 

inertia  about  the  x-axis  when  there  is  no  deformation  of  the  rocket.  In  Eq.  (41), 
I is  a 6x6  matrix,  the  nonzero  elements  of  which  are  1^  * A2+B2  + a^*’i+r2^2* 

*13  " B2  + °r2^1+r2)’  Z22  " *11  * I42  " *24  “ I31  " *13*  I33  “ *44  " B2  + ar22 

T 

and  1^  * ■ 1.  Furthermore,  x ■ ^3  <*>2  “ 3 9 2 9 3^  an<*  C^e  nonzero  elements 

of  the  6x6  matrix  B are 


B12  “ “l[A2  + B2  + ®()l1+r2>2  - W » 


B14  “ fllf2B2'Bl  + 2 a r2(fcl+r2)1, 

»U  * “r2FT  + ai2[B2‘6l  + “VW1 

B16  ‘ VB2  ’ B1  + <’r2(tl+r2)1 


B21  " B12’ 


B23  ’ ‘ B14’ 


B25  ®16’ 


B26  ' B15  ’ 

B32  " *VB2"B1  + 


B33  - -c, 

B34  - 01[2  B2-Bl  + 2or22]  , 

B35  - ur2FT  + ni2^B2"Bl  + or2^  “k* 

»36  ’ “l'VBl  + nt* 


- b32. 


-B34’ 


B44  " “ C’ 


B45  * ‘ B36’ 


(42  cont . ) 


B * B 
46  35 


B53  * B64  " 1* 


where  c and  k are  damping  and  stiffness  coefficients  which  appear  In  the  Internal 


torque. 


The  term  T In  Eq.  (41)  Is  due  to  the  mechanical  thrust  misalignment  find  has 


the  form 


T - [-a  F_  l -a  Fm  l 
y T c z T c 


o]T  , 


where  + (^  + ^ + r2)m2]/M. 


The  complete  set  of  linearized,  rotational  motion  equations  consists  of 


Eqs.  (35a),  (35b),  (36),  (40)  and  (41).  These  equations,  although  linear,  have 


variable  coefficients  if  is  not  constant.  In  the  next  section,  the  additional 


assumption  that  is  constant  is  made  and  closed-form  solutions  for  Qj. 
o>2 > 02»  03»  9 and  ¥ are  obtained. 


2.4  Solution  to  the  Linear  Equations  when  Spin  Rate  is  Constant  and  Nonzero 


spin  i ate,  $ ■ is  constant,  the  matrix  B of  Eq.  (41)  is  constant. 


In  this  case,  the  linear  equations  are  time  invariant  and  may  be  solved  exactly 


by  elementary  methods.  Once  the  solution  for  x(t)  is  obtained,  the  elements 


x^  ■ «2  and  x2  ■ may  be  inserted  into  Eqs.  (35a)  and  (35b)  and  solutions  for 


¥ and  3 obtained  by  directly  integrating  with  respect  to  time. 


I I 
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Me  first  consider  the  homogeneous  equation  lx  = Bx  and  let 


a ■ r1 5 . 


when  is  constant.  Hence,  the  system  of  equations  to  be  solved  is 

x * A x , (44) 


where  A is  a constant  6x6  matrix.  It  is  well  known  that  the  general  solution  to 
Eq.  (44)  is  composed  of  six  linearly  independent  solutions  of  the  form. 


X.  t 


(45) 


where  e^  is  an  eigenvector  of  A and  X^  is  the  corresponding  eigenvalue.  The 


eigenvalues  are  zeros  of  the  function. 


f (X)  - | XI 


(46) 


where  I is  the  6x6  identity  matrix  and 


denotes  the  determinant.  The  eigen- 


vectors e^  are  such  that 


g(Xj)ej  - 0 (J**l,2,3, ...  ,6)  , (47) 

where  B( ) * X^  I - A . 

For  the  particular  problem  of  Interest,  when  j1  0,  the  eigenvalues  and 
eigenvectors  generally  occur  in  complex  conjugate  pairs.  From  this  fact  we 
conclude  that  the  homogeneous  solution  for  each  element  of  x is  the  super- 
position of  six  trigonometric  terms,  each  of  which  is  multiplied  by  an  exponential 
function  of  time. 

We  define  E to  be  a 6x6  matrix,  the  odd  columns  of  which  are  composed  the 
real  parts  of  eigenvectors  corresponding  to  three  eigenvalues,  none  of  which  are 
complex  conjugates  of  another  and  each  of  which  has  a positive  imaginary  part. 


24 


The  even  columns  of  E are  composed  of  the  negatives  of  the  imaginary  parts  of 
the  chosen  eigenvectors.  We  also  let  n.  denote  the  real  part  of  A , and  to 

" J J 

(not  to  be  confused  with  or  <*>3  of  the  previous  subsections)  the  Imaginary  part 
of  Aj.  The  homogeneous  solution  for  x may  then  be  written  as 


-1 


Xh<c)  - E a 5 5o  • 


(48) 


where  x " x(t  ) and,  with  r - t-t  , 
-00  o 


D - 


In  Eq.  (49), 


2n  2 2 

2 222  2 


■33 


Y Y 

e J cos  (OjT  e J sin  to^t 


Y Y 

-e  J sin  WjT  e J cos  to^r 


Equation  (48)  can  also  be  written  as 


(49) 


(j-1,2,3) 


(50) 


5h(t> 


- *(t,t  )x  , 
■ 0-0 


where  t is  the  state  transition  matrix. 

m 

Solution  for  Y and  Q 

To  integrate  Eqs.  (35a)  and  (35b),  we  need  to  replace  0^  and  0^  by  their 
explicit  solutions,  $ by  and  integrate  each  equation  term  by  term.  Dropping 

the  subscript  h for  the  present,  we  write  the  solutions  for  the  x^  in  the  explicit 
forms, 

3 nr 

Xj  - I e 1 (A^  cos  UjT  + sin  u^t)  (j-1,2, . . ,6)  , 


(51) 


where,  if  E.,  is  an  element  of  E and  C is  an  element  of  C * E x , 
J K " J ~ m ~0 


Aij  " E j , 2 i- i °2i-l  + Ej,2i  C2i 


(52a) 


Bij  " Ej ,2i-l  C2i  " Ej ,2i  C2i-1 


(52b) 


Since  ^ = xi  an<*  ^3  - x2»  tlie  desired  expressions  for  these  variables  are 

available  and,  upon  inserting  them  into  Eqs.  (35a)  and  (35b),  we  obtain 

(assuming  $ ■ 4>(t  ) * 0) 

o o 

® * 2 Ji e 1 l(AirBi2)<:os  V 

+ (Ail+Bi2)  cos  q^T 


- (Ai2+Bii)  3in  piT 

- (A^-8!!^  sin 


’ ’ 2 6 1 <<Ai2+Bll)  C°5  P1  + (Ai2-Bll)C0S  V 

+ (A^-8^)  sin  + ^Aii+Bi2)sin 


(53a) 


(53b) 


where  p^  ■ and  qi  ■ Equations  (53)  may  also  be  written  in  the 

matrix  form. 


cos  p t -sin  p.xl  | A . ,-B 


sin  p^i  cos  p^  TJ  LAi2+Bil 


il  12  n^T 
e 


3 rcos  qiT  -sin  qiT1  rA11+Bi2-| 


3 ^ ^ 

+ 7 I 


sin  qiT  cos  q^j  [A^-B^ 


from  which  the  orthogonality  of  0 and  ¥ may  be  seen. 

The  integration  of  Eq.  (54)  is  straightforward  and  the  result  we  obtain  is 


. 3 V 

-7  l ~f~T  l ^ 

1 i*l  n^  + pj 

1 3 .V 
+ i i ~f~r  ^ 

i-l  nj+q2 


A -B 
il  12 


A12+Bil 

Ail+Bi2 


Ai2"Bil 


where 


ni(cpiT_1)  + PiSpiT  -tn1sp1T-pi(cp  t-1)] 
nifspiT*Pi(cpiT_1)1  n1(cpiT-l)+pisp1x 


and  2(x)  has  the  same  form  as  P(r),  but  in  ^(t)  the  qt  take  the  place  of  the  p . 

Equations  (51)  and  (55)  along  with  4>  constitute  the  complete  solution 

for  the  rotation  motion  of  the  two  parts  of  the  rocket  model  when  there  is 
no  mechanical  thrust  misalignment.  To  use  the  solutions  it  is  only  necessary  to 
(1)  determine  the  eigenvalues  and  eigenvectors  of  the  matrix  A,  (2)  form  E, 
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-1 


(3)  invert  E,  (4)  form  c - E x , (5)  determine  A.,  and  B,  .,  (6)  determine  pJ 
■ * -o  ij  ij  i 

and  q^,  (7)  substitute  the  time  at  which  the  solution  is  to  be  evaluated  into 
the  equations,  and  (8)  calculate  the  values  of  the  variables. 

Thrust  Misalignment  Effects 

Because  the  solution  to  the  linearized  equations  for  the  case  when  the 
thrust  vector  is  not  aligned  with  the  x^-axis  will  be  needed  in  the  next  sub- 
section, we  consider  here  the  problem  of  obtaining  such  a solution.  This 
requires  that  we  find  the  general  solution  to  the  nonhomo geneous  equation  (41) , 
or  equivalently,  the  equation. 


where 


x ■ Ax  + f , 


f - I"1  [-a  F_  i , -a  F l .0,0,0.  0]T. 

— » v T c ’ t T n ' 


(57) 


*y  *T  ~c  ’ ~z  *T  ~c  ’ " ’ w ’ w - (58) 

The  solution  for  x(t)  in  this  case  is  composed  of  the  solution  to 
x ” Ax,  l.e.,  the  homogeneous  solution  x^,  plus  a particular  solution.  Since 


f is  constant,  a particular  solution  is 
x “ - A ^f  . 

-p  « - 


(59) 


The  solution  x ■ x,  + x is  subject  to  the  initial  conditions  x * x(0).  Thus, 

- -h  -p  -o  - 

x^(0)  ■ ^0~5p  an8  the  constants  which  appear  in  the  homogeneous  part  of  the  solution 

are  determined  from  Eq.  (48)  with  x^(0)  in  place  of  x^. 

When  thrust  misalignment  effects  are  present,  the  solutions  for  ft^Sx^  and 

Q =x0  have  the  extra  terms  x,  and  x„  and  those  for  0 and  Y have  the  additional 
32  lp  2p 

terms  (x,  /ft.)  sin  +(x.  /ft.) (cos*  -1)  and  (x  /ft  )( 1-cos*)  +(x  /ft  ) sin  *, 
lp  1 Zp  1 lp  1 2p  1 

respectively.  The  effects  of  these  terms  will  be  discussed  in  the  next  subsection, 
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but,  generally  speaking,  their  Inclusion  results  In  periodic  variations  In 
0 and  ¥ at  the  spin  frequency.  Furthermore,  the  inclusion  of  the  effects  of 
thrust  misalignment  in  this  form  (i.e.,  constant  thrust)  results  in  changes 
in  the  bending  amplitudes  and  rates  because  x^(0)  not  x^  is  used  to  obtain 
these. 

For  all  cases  considered  to  date,  one  of  the  complex  pairs  of  eigenvalues 
has  a very  small  magnitude  real  part  and  an  imaginary  part  very  nearly  equal 
to  the  spin  frequency.  Numerically,  this  small  magnitude  real  part  is  sometimes 
negative.  The  other  two  are  always  negative  if  damping  is  present.  That  the 
small  magnitude  real  part  was  sometimes  found  to  be  positive  was  at  first 
surprising,  because  it  implies  that  the  motion  is  asymptotically  stable, 
although  the  spin  is  nominally  about  the  axis  of  minimum  moment  of  inertia.  At 
first  sight,  this  result  appears  inconsistent  with  the  popular  "maximum  moment 
of  inertia  spin  criterion"  of  spacecraft  attitude  dynamics,  which  states  that 
torque-free  spin  of  a quasi-rigid+  body  about  its  axis  of  moment  of  inertia  is 
unstable  when  there  is  inertial  dissipation.11  The  key  difference  in  the  space- 
craft case  and  the  two-body  rocket  is  that  the  spacecraft  must  be  quasi-rigid 
for  the  criterion  to  apply  to  the  two-body  rocket  is  not  quasi  rigid,  since  the 
parts  which  move  relative  to  each  other  are  of  essentially  equal  mass.  To 
verify  that  this  is,  in  fact,  a reason  for  the  different  results,  the  second 
body  of  the  two-body  rocket  model  was  replaced  by  a small  point  mass  attached 
to  a massless  rigid  rod  hinged  to  Body  1 at  point  R.  The  eigenvalues  for  this 
"spacecraft"  model  were  calculated  and  one  of  the  three  complex  conjugate  pairs 
possesses  a positive  real  part.  Hence,  it  appears  that  the  relative  size  of 
the  two  bodies  and/or  lack  of  sufficient  numerical  precision  causes  this  anomaly. 

+The  moving  parts  of  the  spacecraft  have  small  mass  compared  to  the  total 
spacecraft  mass. 
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Actually,  for  our  purposes,  it  is  uot  very  important  if  the  small  magni- 
tude real  part  is  positive,  or  negative,  as  long  as  its  magnitude  is  small 
compared  to  those  of  the  other  real  parts  and  small  compared  to  the  reciprocal 
of  the  time  period  of  interest.  For,  if  this  is  the  case,  the  other  two 
modes  will  rapidly  die  out  leaving  the  lightly  damped  mode  as  the  dominant 
motion. 

The  lightly  damped  mode  is  of  special  interest  because,  generally,  it  is 

essentially  a "rigid-body"  mode  in  which  very  little  relative  motion  between 

the  two  parts  of  the  rocket  exists.  If  the  two  parts  of  the  rocket  are 

identical,  it  is  exactly  a rigid-body  mode.  Furthermore,  the  frequency  of  the 

lightly  damped  mode  is  essentially  the  rigid-body  nutation  frequency  so  that 
• • 

when  0 and  f are  formed  their  parts  due  to  their  mode  are  long-period,  lightly- 
damped  oscillations.  In  a plot  covering  a time  period  much  smaller  than  the 
long-period  of  these  oscillations,  they  appear  as  straight  lines.  Hence,  we 
are  led  to  define  a "steady-state  transverse  rate"  0gs  as  the  square  root  of 
the  sum  of  the  squares  of  the  coefficients  of  the  sine  and  cosine  terms  of 
lightly  damped  parts  of  * x^*  T^ie  square  root  °f  sum  of  the  squares  of 
the  corresponding  coefficients  in  the  “ X2  solution  is  the  same,  i.e., 
and  are  orthogonal. 

The  quantity  0 is  felt  to  be  a good  measure  of  the  adverse  transverse 
ss 

rate  caused  by  bending  and  thrust  misalignment,  the  effects  of  which  are 
incorporated  in  it  through  x^(0)  • xq  - x . 

2.5  Results  from  the  Solution  to  the  Linear  Equations 

Although  the  solution  given  in  subsection  2.4  is  in  closed  form,  it  is 
not  explicit  enough  that  the  effects  of  changes  of  system  parameters  can  be 
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determined  without  numerical  computations.  Hence,  a digital  computer  code 
(FORTRAN  IV)  for  evaluating  the  solution  was  developed.  It  is  completely 

r 

"self-contained"  in  that  no  system  subroutines  are  needed  (except  for  plotting). 

For  the  most  part,  it  is  a linear  system  analysis  code  which  takes  given 
values  of  the  system  constants  and  (1)  generates  the  matrices  I,  B and  T, 

(2)  forms  A * I-1Band  f » A-1  T,  (3)  determines  the  eigenvalues  and  eigenvectors 
of  A,  and  (A)  prints  out  the  explicit  numerical  form  of  the  solution  for  x(t) 
and  those  for  0(t) , V(t),  0N(t)  an(*  ^(t)  at  various  time  points  and  (6)  stores 
these  values  for  plotting.  A working  version  of  this  code  has  been  made 
available  to  U.S.  Army  personnel.  ' 

The  code  was  used  to  obtain  the  results  presented  in  this  subsection. 

These  results  include  "general"  parametric  results  and  some  results  for  a 
specific  rocket  model  which  are  compared  with  test  data.  The  "general" 
parametric  results  were  obtained  by  adopting  a baseline  rocket  model, determining 
its  nominal  behavior  and  then  varying  one  of  the  system  parameters.  The 
- characteristics  of  the  baseline, or  nominal  (N), rocket  model  are  given  in 
Table  1.  The  rocket  model's  two  parts  were  assumed  to  be  homogeneous,  right 
circular  cylinders  of  equal  radius  and  their  moments  of  inertia  were  computed 
from  standard  formulas.  This  method  of  computing  the  moments  of  inertia 
results  values  for  the  overall  moments  of  inertias  which  are  approximately 
the  same  as  those  of  an  actual  free-f light  rocket  . The  length  £j+r^  was  chosen  so 
that  the  center  of  mass  of  the  rocket  model  is  at  the  hinge  point.  This  was 
done  so  the  first  bending  mode  should  be  well  approximated  when  the  two  bodies 
are  identical.  It  is  not  a necessary  requirement. 

■ 
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Table  1.  Physical  characteristics  of  the  baseline  mo^el. 


Body  1 Body  2 


ml 

50.466  kg 

m2 

64.286  kg 

ll 

0.938  m 

*2 

-0.335  m 

rl 

0.938  m 

r2 

0.739  m 

Radius 

0.095  m 

Radius 

0.095  m 

A1 

0.228  kg-m2 

B1 

0.290  kg-m 

*2 

14.91  kg-m2 

B2 

11.86  kg-m' 

44482  Nt 
T 

1T  106.29  kg-m2 
I 0.518  kg-m2 

X 


I 
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The  Initial  conditions  used  to  obtain  all  the  results  given  in  this  sub- 
section were  determined  by  specifying  the  initial  transverse  angular  deflections 
and  rates  of  the  nose  of  the  rocket  and  requiring  that  the  line  segment 
PQ  be  nonrotating  initially  as  it  would  be  in  a launch  from  an  ideally  rigid 
launcher.  Also,  a zero  initial  roll  angle  is  assumed.  These  conditions  are 
sufficient  to  obtain  the  initial  conditions  x^q.  Explicitly, 


and 


’'lO--1(t2+r2)/<tl  + rl)ltVl  + V 1 

x20  * "^*2  + r2^^i  + rl^~0NO^l  + ^NO^  ’ 
x30  " -[(d/U2  + r2)]  X10  ’ 
x40  * “t(<i/(«.2  + r2^  x20  ’ 

X50  " ld,Ul  + V1  ®N0  ’ 

x60  * + V1  h'N0  ’ 

’o  ' “ta2  + r2)/dl  X50 
°o  = -[U2  + r2)/dl  X60  ’ 


(60) 


where  d = + r^  + i + r^ . 

The  initial  conditions  for  the  nominal  case  were  obtained  using  the  above 

equations  and  V„n  » 0.0001  rad.,  0 - i - 0 - 0.  The  time  history  of  f 

NO  NO  NO  NO  N 

is  shown  in  Fig.  6.  The  mean  value  of  V is  very  slightly  decreasing  and 
the  damping  is  apparent  in  that  the  amplitude  of  the  bending  oscillations  is 
noticeably  decreasing  with  time.  The  primary  effect  of  non-zero  initially 
is  seen  in  the  plot  of  ©N  as  a function  of  time  in  Fig.  7.  It  is  a secular  change 
in  of  -1.213x10“ 3 rad/sec.  There  are  no  bending  oscillations  apparent  in 
the  0jj  curve  until  the  Coriolis  coupling  terms  transfer  some  of  the  bending 
present  in  ^ into  the  pitch  plane. 


PITCH  ANCLE  OF  NOSE  (milliradians) 


r 
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TIME  (seconds) 


Figure  7.  Time  history  of  pitch  angle  of  rocket's  nose. 


of  0 were  determined  for  various  spin  rates  In  addition  to  those  listed 


in  Table  2.  For  the  case  of  only  bending  of  magnitude  0.0001  radians,  the  results 


are  shown  in  Fig.  8.  The  variation  of  0 as  a function  of  spin  rate  is 


linear.  When  "additive"  mechanical  thrust  misalignment  (Case  1)  is  also 


present,  the  variation  of  0 with  spin  rate  has  a minimum  at  about  9 Hz 


Optical  lever  data  are  available  for  several  rockets  flight  tested 


Some  of  these  rockets  apparently  bent 


during  an  experimental  program 


0.0001  rad.  Hence,  it  was  decided  to  test  the  adequacy  of  the  two-body 


rocket  model  by  trying  to  match  the  optical  lever  data  of  one  rocket  with 


0 and  ¥ time  histories  generated  using  the  solution 


The  rocket  chosen  to  be  modeled  was  the  GEM  #7.  That  rocket  was  spun  up 


using  an  eroding  spin  turbine,  so  its  spin  rate  was  not  constant.  A mean  rate 


of  47  rad/sec  was  therefore  chosen  because  the  two-body  solution  requires  a 


constant  spin  rate.  The  other  characteristics  of  GEM  #7  are  given  in  Table  3 


The  optical  lever  data  obtained  during  the  flight  of  GEM  #7  were  used  to 


determine  the  "mean"  rates  of  the  yaw  and  pitch  angles  of  that  rocket's  nose 


data,  the  values  of  H*  , 0M,  and  0 due  to  bending  at  the  end  of  guidance 


were  estimated.  These  values  are  as  follows 


0.0  rad /sec  and 


0.00023  radians,  0, 


0.00025  radians,  f, 


0.364  rad /sec. 
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Table  3.  Characteristics  of  GEM  #7 


Length  (m) 

3.348 

Radius  at  aft 

end  (m) 

0.1575 

Mass  (kg) 

114.76 

IT  (kg-m2) 

99.705 

Ix  (kg-m2) 

0.40636 

Thrust  (Nt) 

44482 

Frequency  of 
bending  mode 

first 

(Hz) 

67 

I 


I 
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The  model  characteristics  given  in  Table  4 were  developed  from  the 

actual  characteristics  of  the  rocket.  Note  that  I_  and  I for  the  model 

T x 

are  not  exactly  those  of  the  actual  rocket  because  of  the  use  of  right 
circular  cylinders  to  model  the  two  parts  of  the  rocket. 

The  three  frequencies  of  the  two-body  model  are  82.8  Hz,  68  Hz,  and 
7.44  Hz.  The  characteristics  of  the  solution  are  such  that  the  bending 
frequency  of  the  model  is  68  Hz  and  the  lowest  frequency  is  the  model's 
nutation  frequency. 

To  determine  the  amount  of  thrust  misalignment  (an  equivalent  amount  of 

mass  imbalance  could  be  substituted)  which  results  in  the  mean  rate  cited 

• • 

above,  the  two-body  computer  code  was  run  six  times  with  0NQ,  0Nq, 

and  a ^ successively  equal  numerically  to  0.001.  From  the  six  runs, 

"sensitivity"  coefficients  were  determined.  These,  when  multiplied  by  the 

given  values  of  the  initial  values  of  yaw  and  pitch  angles  and  their  time 

rates  of  change,  and  the  appropriate  thrust  misalignment  angles,  result  in 

the  desired  mean  rates.  The  values  of  a and  a thus  obtained  are 

y z 

0.0006765  and  -0.004598  radians,  respectively.  These  values  are  within  the 
realm  of  physical  possibility. 

By  running  the  two-body  model  computer  program  with  the  values  given 
above,  plots  of  ©N  and  t versus  time  were  generated.  These  plots  were  time- 
shifted  to  match  the  phase  of  the  experimental  data.  A shifting  was  also 
made  to  eliminate  the  bias  of  the  data.  The  resulting  theoretical  time 
histories  are  shown  as  the  solid  curves  in  Fig.  10.  The  dashed  curves  were 
obtained  from  Ref.  12  and  represent  smoothed  optical  lever  data.  The  time 
0.15  seconds  on  the  plots  corresponds  to  the  end  of  guidance  time  of 


PITCH  ANCLE  (MILLTRADIANS)  f , YAW  ANGLE  (MILLIRADTANS) 


Figure  10.  Comparison  of  flight  test  data  and  theoretical  results 
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0.1901  seconds.  Very  good  agreement  Is  apparent,  although  the  model  Is  a 
simple  one.  In  particular,  only  the  first  bending  mode  Is  modeled  and, 
since  a constant  spin  rate  is  used  In  generating  the  analytical  solution, 
the  variable  spin  rate  of  the  actual  rocket  Is  not  precisely  modeled. 

We  conclude  that  the  simple  two-body  model  can  be  used  to  analyze  optical 
lever  data  which  contains  evidence  of  rocket  transverse  bending. 

An  Important  observation  which  may  be  made  concerning  these  results  Is 
that,  on  the  basis  of  the  model,  the  bending  of  the  rocket  contributed  a 
transverse  rigid  body  rate  at  launch  of  magnitude  0.115  rad/sec.  The 
thrust  misalignment  was,  in  this  particular  case,  such  that  It  opposed  the 
effect  of  bending  so  that  the  observed  transverse  rate  at  launch  was  less 
than  that  due  solely  to  bending. 

Table  4.  Characteristics  of  two-body  model  of  GEM  #7 


Body  1 

Body  2 

(kg) 

50.466 

®2 

64.286 

<m) 

0.9376 

*2  (®) 

-0.3354 

^ (m) 

0.9376 

r2  (m) 

0.7393 

Radius  (m) 

0.095 

Radius  (m) 

0.095 

(kg-m2) 

0.2279 

Bx  (kg-m2) 

0.2903 

A2  (kg-m2) 

14.909 

B2  (kg-m2) 

11.857 

k (Nt-m) 

1.491xl06 

Ft  (Nt) 

44482 

1T  (kg-m2) 

106.29 

I (kg-m2) 

X 

0.5182 

42 


2.6  Nonlinear  Equations  for  Motion  on  a Rigid  Launcher 

As  stated  in  Section  1,  the  equations  for  motion  on  a rigid  launcher 
were  not  developed  until  near  the  end  of  this  study.  Because  of  this,  they 
were  not  used  to  obtain  initial  conditions  for  the  free-f light  phase  equations  given 
supra.  This  is  not  considered  to  detract  from  the  results  given  in  the  pre- 
ceding section,  however,  because  the  initial  conditions  used  to  obtain  those 
results  are  representative  of  those  measured  during  flight  tests. 

The  launch-phase  equations  are  given  here  for  the  purposes  of  (1)  "com- 


pleting" the  model  and  (2)  providing  a means  of  studying  the  motion  of  a two- 
body  rocket  ideally  constrained  on  a rigid  launcher.  An  obvious,  and  probably 
needed,  extension  of  the  equations  given  in  this  subsection  would  be  to  include 
a launcher  model  with  three  or  more  degrees  of  freedom. 

Constraints 

While  on  the  launcher,  the  points  P and  Q are  constrained  to  travel  along 
the  x^-axis  of  the  launcher  coordinate  frame 

two  bodies  therefore  have  four  degrees  of  freedom  during  the  launch 
phase.  These  may  be  chosen  to  be:(l)  the  displacement  of  the  point  P along 
the  x^-axis 

(3)  the  pitch  angle  of  Body  1 relative  to  the  Ox^y^z^  frame,  0;  (4)  the  roll 


, Xpj  (2)  the  yaw  angle  of  Body  1 relative  to  the  x^-axis , Y ; 


Ox^yLzL  shown  in  Fig.  11.  The 


angle  of  Body  1 about  the  x^-axis. 

The  orientations  of  the  unit  vector  triads  (iL,  k^) , (i^,  j^,  k^)  and 
(i^,  ^2*  ^2^  re^ative  to  each  other  and  to  the  inertial  reference  frame  unit 


vector  triad  (I,  J,  K)  are  defined  by  the  following  equations: 


where  q is  the  launcher  elevation  angle. 

L 

The  requirement  that  P and  Q translate  along  the  x^-axis  is  fulfilled  if 
P translates  along  that  axis. 


(Il  + -1  + Ei  + *2^  * 1l  * 0 (64a) 

and 

+ £i  + E2  + -2}  * bL  " 0 • (64b) 

Now,  assuming  that  r.  ■ (r  0 0)T,  £.  ■ (X,  0 0)T,  r~  * (r  0 0)T  and 

T I I 1 1 l 

-2  " ^2  0 ^ and  that  02»  ®3*  ^ 311(1  0 are  s™3!1  angles  (which  certainly 

should  be  the  case),  Eqs.  (61),  (62)  and  (64)  can  be  used  to  obtain  the 
following  result: 


(65) 


where  R - (r^  + + r^  + (r2  + By  differentiating  Eq.  (65)  with 

respect  to  time  and  using  the  approximate  equations, 

. (66a) 

* 9c$  + (66b) 

and 

J*3  ■ -6s4>  + ijicij), 

+If  < 0,  then  r„  must  not  be  equal  to  U„  I . 


(66c) 


r 
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we  obtain 


" • “ 

0 

" 0 

"*£2_  " 

2 

. Q 

3 

-R 

2 

• 

9. 

1 

-0. 

£2. 

L 3-1 

L ”2-1 

L 3J 

(67) 


The  total  angular  velocity  of  Body  2 in  the  Body  2 basis  is  denoted  w.  Since 

• • T ^ 

w * A £5  + W9/1  an<*  u2/l  = ®o  S-J  > we  find  that 

2 = I 2 2 


w - (£2, 


(1-R)£22  (1-R)£23)  • 


(68) 


Equations  (65)  and  (68)  are  the  main  constraint  equations  used  in  the 
derivation  of  equations  of  motion  which  follows.  Although  they  are  approximate, 
the  small  angle  approximations  on  which  they  are  based  are  very  good  approxi- 
mations. 

Equations  of  Motion 

Free-body  diagrams  of  Bodies  1 and  2 are  shown  in  Fig.  12.  The  forces  Fp 

and  Fq  are  forces  of  constraint  exerted  by  the  launcher  on  the  rocket,  while 

F_  is  the  force  exerted  on  Body  1 by  Body  2.  By  letting  H . and  m R denote 
-R  ~^,Cj  J 

the  rotational  angular  momentum  and  linear  momentum,  respectively,  of  Body  j 


and  referring  to  Fig.  12,  we  may  write 


Sl/^  * -1  X -R  ‘ -1  X -P  + -1/2  + -s  * 


52/C2  " -2  X -Q  + -2  X -R  ‘ -1/2  ’ 


m,  R„  ■ F_  + m,  g + F + F 
1 — Cp  — T 1 — — p — R 


and 


-2  ?C2  " Zq  “ !r  + m2- 


(69a) 

(69b) 

(69c) 

(69d) 


■ 
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h/h  ■ 'W 


F - r 
-R  -1 


X °1  -c  + £1  X “l  I 
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Also,  by  adding  the  last  two  of  Eqs.  (69),  we  obtain 

M * ?T  + M § + !P  + !q  * <72> 

where  as  before  and  R^  is  the  acceleration  of  the  system  center  of  mass. 

If  we  neglect  the  frictional  forces  which  act  parallel  to  the  x^-axis  at  P and 
Q,  then 

M \ m h ■ k * H ? *L  (73) 

c 

AAA 

Letting  F * F i + F j + F k and  using  the  small  angle  approxima- 
tion  to  A to  get 


-2  + -2  “ ^2+r2^  ^il  + 0i  li  " e2  -l^  * 

mi  may  rnnhinp  F.qa,  fTH  tO  get 


h/CL  2 I1/(R-1)]52/C2  + <W\<*2  ll  + e3  *1> 

-[1/(R-1)][12  x(m2  R - m2  g) ] 

- r^  Rc  - mx  g]  + [R/(R-l)]T1/2  + Tg  (75) 

At  this  point,  we  begin  the  change  to  matrix  notation  and  express  R_  and 

+ i 

R in  their  Body  1- basis  forms, 

"C2 


The  brackets  with  subscripts  indicate  that  the  inertial,  or  total,  vector 
derivative  bracketed  is  evaluated  in  the  basis  corresponding  to  the  subscript. 
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and 


[R„  ] - B(x  0 0)T  + nSr,  -r.il 

"C1  1 " P III1  I1  I 

[rc  )1  - B(kp  o 0)T  + 


(76a) 


.T  -- 


.T  - 


+ A uui  r„  - A r2  w 
22  2 "22 


Nov,  from  Newton's  second  lav,  still  using  vectors, 


Rc  - !t/m  + g 

and, from  kinematics  and  the  definition  of  the  center  of  mass. 


(76b) 


(77) 


M Rc  " M<5p  + £i>  + m2(Ii  + (78) 

Thus,  by  expressing  the  derivatives  of  the  vectors  in  Eq.  (78)  in  the  launcher 
basis  using  Eqs.  (76),  we  get 

(x_,  0 0)T  - E BT  {Ft/M  - nCr,  + r,  Q 

-P  L -1  - jT  ml  i1  I 


-g[(0fl  4,  - i,  0)  + A (bKd  r,  + (R-l)r,  n] 

II  I II  22  " 2 I 

+ C(0  0 g)T}, 


where  we  have  used  the  constraint  equation  (68)  to  write  u 

• T T 2 

(l-R)Qj)  and  r ■ (r2  0 0)  has  been  assumed,  also 


(79) 


( Qj_  (1-R)  02 


§1 


1 

0 

L0 


0 

0 

0 


0 
0 
0 J 


49 


We  next  use  Eqs.  (76)  and  Eq.  (79)  in  Eq.  (78),  expressed  in  matrix  form, 
to  obtain  the  equation  which  governs  the  rotational  motion  of  Body  1, 

LnN  a - - Q I-  Q + AT  w I-  to [ 1/  (R-l)  ] 

ON  1 X 1^"  1 2 2^  2 

-p(il1+r1)(FT  - M g sin  0L)  (0  02  02)T 

- m9[l/(R-l)]AT  1 (A[B  E BT(F_/m  - 00  r 

i 1 II  I 

- vi(£KU,  + cow  r_)  + C(0  0 g)^ 

III  22  2 

+ 00(£..+r, )]  + ww  r»  - A C(0  0 g)^} 

” -1  -1  22  22 

- m,  r B E BT  [FT/M  - 00  r. 

TT  jj  X1 

- p(00  l.  + AT  W(L  r2)  + C(0  0 g)T] 

II  I 22  " 

- m1 [r.  00  r - r.  C (0  0 g)T] 

1 I II  I I 

+ [R/ (R-l)  ]A(0  + c @2  k 03  + c 03)T 

+ (T  0 0)T  , (80) 

where 

■ON  il  + ■ *2  »-mi^l-l  “ £l  2 Si  2 

- [1/  (R-l)  ] m^tA1  \2  AU^)] 

- &T  l2  2 2 2i  2TCl1  + - AT  r2[l/(R-l)]) 

-(m^)  rx  B Ex  B1^  - AT  \2  [1/(R-1)  ] J>  (81) 

and  the  axial  component  of  FR,  FR  , is  approximated  by  y(FT  - m 8 sin  ©L) • 


Numerical  integrations  of  Eqs.  (66),  (79)  and  (80)  can  be  accomplished 
via  digital  computer,  using  Eqs.  (65)  and  (67)  to  compute  6^  and  their 
time  derivatives  when  required. 

2.7  Results  for  Motion  on  a Rigid  Launcher 

As  an  example  of  the  results  which  may  be  obtained  by  numerically 
integrating  the  on-launcher  equations  of  motion,  the  solutions  for  ©N  and  4! 
as  functions  of  time  for  the  GEM  #7  model  are  given  in  Figs.  13  and  14.  To 
be  noted  is  that  the  initial  value  of  0N»  which  is  due  to  the  bending  of  the 
rocket  caused  by  its  weight,  is  small  in  magnitude.  As  the  thrust  is 
applied  and  the  rocket  spins  up  (an  eroding  spin  turbine  is  simply  modeled 
by  using  Ts  ■ 0.02  FT) , the  pitch  angle  of  the  nose  oscillates  about  a mean 
value  of  about  8.8xl0~5  radians  with  damping  apparent.  The  transfer  of 
bending  into  the  yaw  plane  through  the  Coriolis  and  Euler  acceleration  terms 
is  also  apparent  in  Fig.  14. 

An  important  fact  that  requires  further  study  is  that  the  bending  due  to 
gravity  and  thrust  is  not  as  large  in  magnitude  as  that  observed  in  certain 
flight  tests;  for  example,  those  of  GEM  // 7. 


TIME  (SECONDS) 


Figure  13.  Theoretical  attitude  of  rocket's  nose  during 
launch  - pitch  angle. 
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Figure  14.  Theoretical  attitude  of  rocket's  nose 
during  launch  - yaw  angle. 
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SECTION  3.  CONTINUOUS  ROCKET  MODEL 
3.1  Description  of  the  Model 

The  continuous  rocket  model  developed  during  this  study  Is  an 
extension  of  that  described  In  Ref.  9.+  The  Ref.  9 model  Is  that  of 
a pinned-plnned-free  slender  beam,  which  has  Its  points  of  support 
constrained  to  move  in  a fixed  straight  line  during  launch  and  free 
to  move  during  free-flight.  The  results  obtained  using  this  pinned- 
plnned-free  model  did  not  agree  very  well  with  optical  lever  data. 

Hence,  the  present  model  was  conceived. 

The  continuous  model  actually  is  two  models  In  the  sense  that  during 
on-the-launcher  motions  of  the  rocket,  its  bending  deformations  are 
described  using  mode  shapes  for  a pinned-plnned-free,  uniform,  slender 
beam,  as  in  the  previous  model.  However,  during  free-flight  motion  of 
the  rocket,  its  bending  is  modeled  using  free-free  mode  shapes  for  a 
uniform,  slender  beam.  The  use  of  these  mode  shapes  in  the  equations  of 
motion  presented  here  iss  not  necessary,  since  mode  shapes  obtained 
experimentally,  or,  for  example,  by  using  finite  element  computer  codes, 
could  be  used.  In  the  absence  of  such  mode  shapes,  however,  these 
approximate  ones  were  used  to  obtain  results  given  later  in  this  section. 

Sketches  of  the  physical  model  of  the  rocket  during  the  launch  phase, 
when  two  points  of  the  rocket  (points  P and  Q)  are  constrained  to  travel  in 
a fixed  straight  line,  and  during  the  free-flight  phase  when  the  rocket's  motion 
is  unconstrained  are  shown  in  Fig.  13.  The  rocket  is  assumed  to  be  a linearly 

number  of  misprints  are  present  in  Ref.  9 . A corrected  copy  can 
be  obtained  from  the  first  author  of  this  report  upon  request. 
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elastic  body  in  the  sense  that  it  may  deform  transversely;  however,  it  is 
assumed  to  be  rigid  longitudinally.  A small  segment  of  the  rocket,  obtained 
by  cutting  it  transversely  at  two  closely  neighboring  longitudinal  stations 
is  considered  to  have  symmetric  elastic  properties,  but  the  stiffness  and 
mass  density  of  the  rocket  may  be  functions  of  longitudinal  position. 


Figure  15.  Sketches  of  continuous  rocket  model, 
(a)  on  launcher,  (b)  off  launcher. 


j 

! 
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The  mass  of  the  rocket  is  assumed  to  be  constant,  although  the  thrust 
force  which  is  assumed  to  act  on  it  is  generated  by  the  mass  flow  of  burned 
propellant  from  it.  This  assumption  appears  justified  if  the  model  is  used 
in  modeling  motion  of  a rocket  during  only,  say,  the  first  second  of  its 
flight,  provided  its  mass  changes  by  only  a small  percentage  during  that  time 
period . 
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Only  the  external  force  due  to  gravity,  an  "external"  force  and  moment 
due  to  thrust  and,  during  the  on-launcher  motion,  constraint  forces  are 
assumed  to  be  present.  The  thrust  force  is  assumed  to  act  at  the  front  end 
of  the  "motor  section"  of  the  rocket  as  shown  in  Fig.  16  in  determining  the 
internal  force  due  to  thrust.  However,  for  the  purpose  of  computing  the 
force  and  torque  on  the  rocket  due  to  thrust  misalignment,  the  thrust  is 
assumed  to  act  on  the  aft  end  of  the  rocket.  These  assumptions  are  compatible, 
since  the  misalignment  torques  and  forces  are  generated  by  undesirable 
turning  of  the  flow  from  the  motor  due  to  the  rocket's  nozzle  not  being 
perfectly  aligned  mechanically  with  the  rocket's  longitudinal  axis  and  to  the 
rocket's  bent  shape,  while  the  principal  axial  force  should  be  considered  to  act 
at  the  front  end  of  the  motor  section. 8 


Figure  16.  Positioning  of  the  thrust  force. 


In  Fig.  17,  three  coordinate  frames  used  in  developing  equations  of 
motion  are  shown.  The  OXYZ  frame  is  fixed  to  the  launcher  and  is  assumed  to 


be  an  inertial  ("rigid"  launcher)  frame.  The  Pxyz  frame  has  its  origin  at  the 
aft  end  of  the  rocket  with  its  x-axis  passing  through  the  point  Q.  It 
rotates  in  roll  with  the  rocket.  The  Cuj^u^  frame  has  its  origin  at  the 
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center  of  mass  of  the  rocket  and  the  u^-axes  are  instantaneous,  centroidal, 
principal  axes  of  the  rocket.  The  Pxyz  system  is  of  principal  Importance 
during  the  launch  phase,  since  then  the  points  P and  Q are  constrained  to 
move  in  a straight  line.  During  the  free-f light  phase,  the  Cu^U2u3  system 
is  used  as  the  reference  frame  with  respect  to  which  bending  deformations 
are  measured  because  the  equations  of  motion  are  thereby  simplified  and  be- 
cause the  free-free  mode  shapes  of  a uniform  beam  due  to  transverse  bending 
may  be  logically  measured  with  respect  to  such  a frame.  Because  different 
coordinate  frames  are  used  for  the  launch  and  free-f light  phases,  coordinate 
transformations  between  the  frames  are  needed.  These  are  discussed  following 
derivations  of  the  required  equations  of  motion  for  each  phase. 


Figure  17.  Coordinate  frames  used  in  analysis  of  continuous  model. 

Ue  should  point  out  that  in  this  section  some  of  the  same  symbols  used  in 
Section  2 are  used  in  referring  to  different  physical  quantities  than  before. 
Hence,  each  section  should  be  considered  separately. 


3.2  On-Launcher  Motion 

For  motion  during  the  launch  phase,  we  need  to  determine  the  motion  of 
the  Pxyz  system  and  the  motion  of  the  rocket  relative  to  that  system.  To  do 
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this,  we  consider  the  arbitrary  rocket  element  of  length  Ax  shown  in  Fig.  18 
In  that  figure,  r(x)  is  the  position  of  the  aft  end  of  the  element  with 
respect  to  point  P,  -f(x)  is  the  intemal+  force  acting  on  the  -x-end  of  the 
rocket  element  and  -m(x)  is  the  moment  on  that  end.  The  vectors  r and  r 

• " E — D 

locate  the  centroids  of  the  -x-  and  +x-ends,  points  a and  b,  respectively. 
The  mass  per  unit  length  of  the  rocket  is  a(x)  and  a(x)gAx  is  the  gravita- 
tional force  on  the  element.  As  shown,  the  force  and  moment  on  the  +x-end 
of  the  element  are  obtained  from  those  on  the  -x-end  by  reversing  their 
senses  appropriately  and  using  the  first  two  terms  of  the  Taylor  series 
expansions  of  the  resulting  force  f(x)  and  moment  m(x)  about  x. 


The  vectors  r and  r.  may  be  approximated  by 
-a  -b 


rft  = -(3r/3x)Ax/2 


and 


rfe  - (3r/3x) Ax/2, 


(82a) 


(82b) 
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respectively.  And,  since  the  moment  about  c is  given  by 
m - r x(-f(x))+  r,  x(f(x)  + 3f(x)/3x  Ax) 

C “3  — D ” “ 

+ 3m(x) / 3x  Ax  , (83) 

through  first  order  in  Ax, 

2?c  = { [(3r(x)/3x)x  f(x)]  + 3m(x)/3x}Ax.  (84) 

The  net  force,  fc>  on  the  element  can  also  be  approximated  by 

f 2 [ (3f/3x)  + a(x)  g]  Ax,  (85) 

By  letting  Ax  -►  dx,  Eqs.  (84)  and  (85)  become  exact  in  the  differential 
sense.  Then,  by  letting  a£  denote  the  acceleration  of  the  center  of  mass  of 
the  rocket  segment  a(x)dx,  we  have 

aac  * 3f/3x  + ag  • (86) 

Also,  if  we  let  h^  denote  the  angular  momentum  of  the  element  about  c, 

hdm  * [(3r(x)/3x)x  f(x)]  + 3m(x)/8x.  (87) 

Now,  we  may  write  r(x)  as 

r(x)  ■ x j + y(x)  J + z (x)  k , (88) 

so  that  , 

3r(x)/3x  - i + 3y/3x  j + 3z/3x  k . 


(89) 
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e„  *-Oc/3x  + 3z  /3x)  , 

L O 


(93b) 


6,  * 3n/3x  + 3y  /3x  , 
J o 


(93c) 


0 - (0  0 0 ) , 
- x y z 


(93d) 


a - (o  e2  e3)  , 

m 

02  - - 32?/3x3t 


(93e) 


(93f  ) 


03  * 32n/3x3t 


(93g) 


It  follows  therefore, from  (vector  notation)  h,  * (6h,  /6t)  + (fi+A)x  h, 

— am  —am  m — — —am 

and  Eq.  (87),  that 

d^(A  SJ  - A £ n + A)  +(A  Q + A)d^(A  0 + A) 

mm  m m 

» + A[(3r(x)/3x)  f(x)  + 3m(x)/3x]  . 

The  acceleration  ac  may  also  be  expressed  in  the  matrix  form, 


a ' B R + 20r  + nnr-r0  + r , 
-C  ■ LP 

where 


B ■ 0 c$  3$ 

.0  — s 4>  c$ 


R is  the  acceleration  of  point  P,  r ■ (0  3n/3t  Sij/at)1  and 

LP  T 

r-(0  32n/3t2  32C/3t2)T  . 


Using  Eqs.  (94)  and  (95)  in  Eq.  (90)  and  neglecting  the  moments  due  to 

■f* 

components  of  f(x)  other  than  the  axial  one,  which  we  now  call  f(x),  we  get 

UTdj(A  Q - H Q + i)  + AT(A  0 + A)dj  (A  0 + X)} 

m m id  m ni 

+ j^(0  02f(x)  93f(x))T  + 32m(x)/3x2 

- (0  aa  -oa  )T  + (0  f - f )T  - 0 , (97) 

c c ft  ft  — 

z y Bz  6y 

where  a and  a are  the  y-  and  z-components,  respectively,  of  the  acceleration 

Cy  Cz 

of  c and  f and  f are  analogous  components  of  og.  The  term  in  curly 
gy  gz  2 

brackets  is  the  time  rate  of  change  of  the  rotational  angular  momentum  of  the 

rocket  element  about  its  center  of  mass.  If  the  rocket  has  an  essentially 

uniform  mass  distribution  and  is  very  slender,  this  term  should  vary  little 

with  position  x.  Hence,  to  simplify  matters  considerably,  we  neglect  the 

partial  derivative  of  the  bracketed  term  with  respect  to  x.  We  then  use  the 

T 

fact  that  on  the  launcher  ft  ■ (ftx  0 0)  and  also  use  Euler's  bending  moment 

equation  to  obtain  the  following  equations  which  govern  the  bending  of  the 
rocket  while  it  is  on  the  launcher: 

of  n-  2«x;  - «x(*0+c)  - «x2(yo  + n)l 

— [ (y * + n’)  f(x)]’  -f  + (E  I n")"  - 0 (98a) 

o gy 

ofii  + 2ftx  n + ftx(yQ  + n)  ~nx2(z0  + 

-f(z*  + c')  f(x)]'  -f  + (E  I ;")"  - 0,  (98b) 

o gz 

+The  force  f(x),  as  defined  in  Appendix  A,  has  a transverse  component  only 
when  x*0. 


I 


where  ( )-3(  )/3t  and  ( )'-3(  )/3x. 


Equations  (98)  are  subject  to  the  following  boundary  conditions: 


At  point  P: 


Zero  Deflection 


n - c - 0 


Zero  Moment 


n"  - C"  - 0 . 

At  the  nose  of  the  rocket,  x - L: 


(2  conditions) 


(2  conditions) 


Zero  Moment 


n"  - c"  - 0 


Zero  Shear 


(E  I n") ' -(El  ?")’  - 0 


(2  conditions) 


(2  conditions) 


Also,  at  the  point  Q,  n and  z,  must  be  zero  (continuous)  and  have  continuous 
first  derivatives  with  respect  to  x.  These  are  four  more  conditions. 

The  equations  which  govern  the  rocket's  center  of  mass  and  its  rotation 
about  its  center  of  mass  while  it  is  on  the  launcher  may  be  obtained  from 
Equations  (86)  and  (87) , by  integrating  those  equations  over  the  length  of 
the  rocket.  From  Eq.  (86),  we  have 

*vt-  L L 

f era  dx  ■ f Sf  W + fsffxWSvdx  + f eto  H v 


Joa  dx  - J W3f(x)/3xdx  + /3f(x)/3xdx  + joa  dx, 

^'C  ° “ *M  ' ° 


where  the  integral  over  3f(x)/3x  has  been  broken  into  two  pieces  to  allow  for 
the  discontinuity  in  f(x)  at  Now.  the  left-hand  side  of  this  equation 
is  the  total  mass  M^  times  the  acceleration  of  the  center  of  mass  a^,.  Also, 
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neglecting  any  frictional  forces  at  P and  Q and  referring  to  the  definitions 
of  f(x)  given  in  Appendix  A,  we  get. 


“l  *C  “ -T  + MT  ^ 


(100) 


We  are  concerned,  during  the  launch  phase,  with  only  motion  along  the 
x^-axis.  Hence,  we  need  only  the  x^-component  of  Eq.  (100).  This  can  be 


written  as 


Hj  v»L  - FT  - Hj,  g sin  @L  , (101) 

P 

where  u^  is  the  x^-component  of  the  velocity  of  P and  0^  is  the  launcher 
P 

elevation  angle. 

While  it  is  on  the  launcher,  the  rotation  of  the  rocket  is  constrained 
in  such  a manner  that  only  rotation  of  the  Pxyz  frame  about  the  x-axls  is 
allowed.  Although  the  bending  of  the  rocket  produces  variations  in  the  moment 
of  inertia  about  the  x-axis  they  are  very  small  compared  to  the  axial  moment 
of  inertia  of  the  undeformed  rocket.  Thus,  letting  T^  denote  the  x-component 
of  the  external  torque  acting  on  the  rocket. 


0 - T /I  . (102) 

X XX 

In  addition  to  the  Eqs.  (101)  and  (102)  the  kinematic  equations, 

*L  " \ (103a) 

P P 

and 

♦ - 0 , (103b) 

x 


are  needed. 
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Equations  (101)  through  (103)  are  uncoupled  from  Eqs.  (98).  The  latter 
are  partial  differential  equations  with  time  varying  coefficients  and 
apparently  can  only  be  solved  exactly  by  using  finite  difference  techniques. 
Alternatively,  approximate  solutions  to  Eqs.  (98)  may  be  obtained  by  obtaining 
a set  of  ordinary  differential  equations  using  the  method  of  assumed  modes.10 
This  alternative  approach  Is  attractive  and  also  produces  accurate  results 
if  the  mode  shapes  are  chosen  carefully  and  an  adequate  number  of  them  is 
used. 

The  mode  shapes  chosen  for  use  herein  s re  those  for  the  transverse  vibration 
of  a pinned-plnned-free,  uniform,  slender  beam.  These  were  also  used 
in  Ref.  9,  but  there  they  were  used  for  both  the  launch  phase  and  free-f light 
phases  of  motion, while  here  they  are  used  only  during  the  launch  phase.  How 
these  mode  shapes  ^ are  determined  is  discussed  in  Appendix  B.  They  are 
comparison  mode  shapes,  in  that  they  satisfy  the  geometric  boundary  conditions 
and  also  the  conditions  of  zero  bending  moment  at  each  end  and  zero  shear  at 
the  nose  of  the  rocket. 

To  transform  Eqs.  (98)  into  ordinary  differential  equations  we  approximate 
n and  4 by 


n “ l P1 
j-1  3 3 


(104a) 


t q*  *4  » 


(104b) 


respectively,  where  the  p^  and  q^  are  functions  of  time  t and  the  ^ are,  of 
course,  functions  of  x only.  By  substituting  Eqs.  (104)  into  Eqs.  (98), 


multiplying  each  of  Eqs.  (98)  by  ip^,  integrating  over  the  length  of  the  rocket 
and  performing  a good  bit  of  algebra,  we  obtain  the  matrix  equations. 


M(p  ~2Qx  q - &xq  - Q*  p)  + (K  - FT  £)p 

-(g  sin  4>)m  + FT£-^xb-n^a  + Dp-0 

M(g  + 2flx  g + fix£  - q)  + (K  - Ft  r)q 

-(g  cos  $)m  tF^d+^a-a^b+Dj'O 


(105a) 


(105b) 


where  M is  an  NxN  matrix  with  elements. 


m±  - / a dx  ; 

J 0 J 


(106a) 


£ and  q are  Nxl  matrices  with  elements  p^  and  q ^ , respectively;  K and  £ are 


NxN  matrices  with  elements, 


k - J MEI  *")"  dx 

ij  0 J 


Yi1  “ J "VV  f(x)  dX  ’ 

J 0 J 

respectively;  m is  an  Nxl  matrix  with  elements, 
L 

m.-J  crip  dx  . 

1 0 1 


(106b) 


(106c) 


(106d) 
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Also,  a,  b,  c and  d are  Nxl  matrices  with  elements, 
L 


and 


ai  " i"  ° *i  y0  dx 
0 

L 

b * / o i|<  z dx 
* e\  1 o 


C,  “ J y'  f(x)  dx 

1 o ° 

L 

di  ■ J ^ 2^  f(x)  dx  , 
0 


(106e) 
(106f ) 

(106g) 

(106h) 


respectively.  Finally  D is  an  NxN  damping  matrix.  Equations  (102) 
through  (103)  and  (105)  are  the  governing  equations  for  the  launch  phase. 

Since  the  principal  type  of  attitude  data  available  for  free  rockets  is 
optical  lever  data  and  since  the  quantities  measured  are  optical  lever  are 
the  pitch  and  yaw  angles  of  the  nose  of  the  rocket,  we  require  expressions  for 
these  angles  in  terms  of  the  roll  angle  $ and  the  p^  and  q^ . Because  the 
body-fixed  angles  the  rocket's  nose  makes  with  the  x-axis  are  6 *-(3?/3x) 

x**L 

and  Hence» 


and 


02N  3 “l 

j-1  J J 


(107a) 


(107b) 


-ar-  * --  - 


These  may  be  transformed  into  pitch  and  yaw  angles  0^  and  4^  of  the  rocket's 
nose  using  the  expressions. 


0 m 0 +9  cos  $ — 9 sin  $ 

N L 2N  3N 


(108a) 


and 


’’h  * e2S  Sl“  ‘ + % C0S 


(108b) 


3.3  Free  Motion 

In  describing  the  free  motion  of  the  rocket,  we  switch  to  the  coordinate 
frame  was  defined  previously.  If  the  rotary  inertia  term 

3 h^jj/3*  is  again  neglected  and  the  internal  force  treated  as  before,  the 
bending  equations  for  the  free-flight  phase  ate 


-(E  I u2fi) " + [*(*)  u2^'  + °§'i 2 + *T  ^3p  + az* 


*§2  + 52e  -2“l“3e  S u2  + Vl"l 


- “l  u3  + “3  "l1 


(109a) 


and 


-[E  I u'^]"  + [f(x)u^]'  + 08*«3  + Ft  (92p  - <0  6(x-0) 
- o[R  -e3  + u3e  + 20,^  - u3  + «3»lUl 


+ a)^  - u>2u1]  , 


(109b) 


where  u2e  and  u^  are  elastic  deformations  in  the  u2~  and  u^-directions , 
e2  and  e^  are  unit  vectors  with  the  u2~  and  u^-axes,  R£  is  the  acceleration 
of  the  center  of  mass  of  the  rocket,  is  the  u^ -component  the  angular 
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velocity  of  the  coordinate  frame  and  Rc  is  the  acceleration  of  the 

center  of  mass  of  the  rocket. 

Nov,  since  only  the  external  forces  of  thrust  and  gravity  are  assumed  to 
act  on  the  rocket. 


Rc  - + g , (110) 

Hence,  the  gravity  terms  disappear  from  Eqs.  (109)  when  Eq.  (110)  is  used  and 
we  get 

o[u2e  -2  oo^u.^  _U)^U2  + U1<*>2U1  “ (*)lu3  + u3ul^  + FT^93p  + “z^  ^(x-0) 

- -[E  I u^]"  + [f(x)u’l’  - (o/MT)FT-e2  (111a) 


and 


a[u3e  + 2u,1u2e  - u2u3  + o)1u3ul  + Ulu2  - u)2u2]  + FT  02p  - ay)  5(x-0) 


--[El  u'^]"  + [f(x)u’]’  - (o/VIt‘^3 


(111b) 


As  stated,  the  axes  of  the  CU^U2U3  sFsCein  are  at  all  times  the  centroidal 
principal  axes  of  the  rocket.  Ue  shall  use  this  fact  now  in  deriving  equations 
which  govern  the  rocket's  rotational  motion. 

The  acceleration  of  the  arbitrary  element  of  mass  dm  of  the  rocket  shown 
in  Fig.  19  is 


Sdm " 5c + “ + £ » U12) 

where  u is,  as  before,  a vector  from  C to  a point  c on  the  elastic  axis  and 
p is  a vector  from  c to  the  mass  element.  The  point  c is  also  assumed  to  be 


the  center  of  mass  of  segment  of  the  rocket  of  length  du^» 
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L?*C  L/XC 

y u X u odux  +y  (dj  • idm+  o;dm  X dj 

-x„  -x„ 


2*»>  - I • 


(115) 


where  xc  is  the  x-coordinate  of  the  center  of  mass  of  the  rocket  in  the  Pxyz 
frame.  Now,  ujdm  is  equal  to  the  angular  velocity  w plus  the  angular  velocity 
of  the  segment  of  the  rocket  about  c.  For  a slender  rocket,  the  contribution 
to  the  left-hand  side  of  Eq.  (115)  due  to  the  angular  velocity  of  the  segment 
about  its  center  of  mass  is  much  smaller  than  that  due  to  the  remaining  terms. 
Furthermore,  to  neglect  this  contribution  is  consistent  with  the  approximation 
of  zero  rotary  inertia  in  the  bending  moment  equations.  Thus,  we  obtain 


• 

I u x u adu^  + I (d^  • u + aj  x dj  • w)  = T 
-x„ 


______  _ (116) 

”C  ‘XC 

Finally,  we  shall  neglect  the  small  terms  in  the  second  integral  in  Eq.  (116) 
which  are  due  to  the  rotation  of  the  segment  of  the  rocket  with  respect  to 
the  Cu^u^u^  frame  due  to  bending.  This  is  also  consistent  with  the  assumption 
of  a slender  rocket. 

Now, 


ti  ■ °u°  + 2o)XU  + uxuixu  + uxu  , 


(117) 


and  because  the  u^  axes  are  principal  axes  we  can  use  results  of  that  fact, 


L-xc 

^"crUjUj  du^  - 0 (i,j-l,2,3,ii*j)  , 

-x„ 


(118a) 


C 

L-x, 


f GUjUj  duL  - 0 (j-2, 

-x„ 


3) 


(118b) 


and 


J 
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L"XC 

J Ou^  duL  « 0 (j-2,3). 


(118c) 


Co  obtain  from  Eq.  (116)  Che  result. 


T*jS+a>xI  *u-T, 


(119) 


where  1 is  the  centroidal,  principal  inertia  dyadic.  The  elements  of  I must, 

of  course,  differ  slightly  from  the  principal  moments  of  inertia  of  a 

perfectly  straight  rocket,  but  such  differences  are  second  order  in  u^q  and/or 

u • Hence,  they  are  neglected.  It  follows  that  the  only  way  in  which 
Je 

Eq.  (119)  can  be  affected  by  the  bending  of  the  rocket  is  through  the  torque 
T. 

The  torque  T is  due  to  misalignment  of  the  centerline  of  the  nozzle  with 

the  the  u^-axis.  As  shown  in  Fig.  20,  this  misalignment  is  assumed  due  to 

three  factors;  the  mechanical  misalignment  of  the  nozzle  (a  and  a ),  the  non- 

y z 

straightness  of  the  rocket  [P ( p and  (^U2Q^ui^p^  and  bending  of 
the  rocket  [(3u3e/3ui^P  and  (3u2e^3ul^‘  The  tor9ue  is 


T - FT  xc  {tOu3/3Ul)p  - oy]  e2 

- [(3u2/3Ul)p  + otz]  e3  j , 


(120) 


+The  u.-  are  the  counterparts  of  y and  z 
JO  o c 


Figure  20.  Orientation  of  thrust  vector. 

The  translation  of  the  center  of  mass  of  the  rocket  is  governed  by 
Eq.  (110);  however,  we  chose  to  rewrite  it  in  the  form. 


V + ui  x V 


F/M 
-T  T 


+ g 


(121) 


where  V is  the  velocity  of  C with  components  u,  v and  w and 

?T  * V^l  + [<3u2/3u1)P  + S2 

-[-(3u3/3u1)p  + oy]  e3]  . ((122) 

Finally,  the  kinematic  equations  for  translational  and  rotational  motion 
are  the  usual  ones. 


e 
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and 


= +(«  sd>  + c$)tan  0 

(124a) 

* u>2  c<{>  - u>3s<|> 

(124b) 

* (to 2 s<j>  + OlJ^ c<p)  / cos  0, 

(124c) 

respectively. 

Equations  (111) , (119) , (121) , (123)  and  (124)  are  those  required  to 
completely  specify  the  bending  of  the  rocket  and  the  motion  of  its  center  of 
mass.  However,  we  still  must  transform  Eqs.  (Ill)  into  ordinary  differential 
equations  if  we  are  to  solve  those  equations  in  the  same  manner  as  the  launch- 
phase  bending  equations. 

As  stated  at  the  first  of  this  section,  we  shall  use  the  mode  shapes  for 
a slender,  uniform  free-free  beam  to  make  the  required  transformation.  The 
mode  shapes  are  denoted  <t>^  and  are  functions  of  x.  However,  it  is 
obvious  that  u^  * x-x^,.  They  are  discussed  in  Appendix  C. 

Proceeding  as  in  the  launch  phase  case,  we  set 


2e 


and 


■ I Vj 

j«i  -1  3 

H 

* l Qj  *1  • 

j-1  3 J 


(124a) 


(124b) 


-^e  Pj  and  are  time  functions.  By  inserting  these  approximations 
(i*  III i,  performing  the  necessary  multiplications  by  ^ and  inte- 


«•  .ength  of  the  rocket,  we  get 
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Mf(P  - 2 ux  Q ■ ixQ  - uj  P)  + (Kf  - Ft  rf)P 

+ (^3  + ” “l  -f  “ wi  -f  + Ft  -f  ■FT^03p+az^^0^ 

+ (F^H^  rnf  - 0 (125a) 

and 

Mf(Q  + 2 Ul  P + (Lx  P - u2  Q)  + (Kf  - Ft  rf)Q 

+<“l“3  " V-f  + “l  -f"ul  -f  ‘ FT  -f  •FT(02p"ay)^(O) 

+ <ft/mt>3  ?f  " 0 * (125b) 

where  the  subscript  f is  used  to  indicate  that  the  matrices  so  subscripted 
are  for  the  free-flight  phase.  The  matrix  is  NxN  with  elements 


*4 

/°  V] 


(126a) 


The  matrices  P and  Q are  Nxl  matrices  with  elements  P^  and  , respectively. 

Also,  the  matrices  K,  and  r,  are  NxN  with  elements, 

■f  -f 


w 

kijf  '/vE  1 ♦?" d* 


(126b) 


’ f*i*l 


f(x)  dx, 


(126c) 


respectively.  Furthermore,  the  elements  of  the  matrices  a^,  b^,  c^,  d^, 
and  m^  are  defined  as  follows: 


‘i.  ' / 0*1  “j 


(126d) 
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>1,  ■ / °*i  “30  d*  ' 


=3  ' / ♦*' 

£ 0 

‘ ■ / *1  “30 


f (x)  dx 


f(x)  dx 


l • / ♦!<*-«<;)  , 


u 

V ■ / V dx  » 


and 


<VV2  ‘ <WIuio  + ^1P3*3(0>  + “.1 

Wl  ■ <Wtu30  + jx  03  »]<0)  - «yl 


(126e) 

(126£) 

(126g) 


(126h) 

(1261) 

(126j ) 

(126k) 


3.4  Transformation  Required  at  Launch 

Because  different  reference  frames  and  mode  shapes  are  used  In  describing 
the  motion  of  the  rocket  during  the  launch  and  free-free  phases,  a rather 
complicated  transformation  Is  required  at  the  Instant  of  launch.  Simply  put, 
we  must  determine  "Initial"  conditions  on  the  free-flight  variables  using  the 
"final"  conditions  on  the  launch  variables.  Here,  of  course,  the  instant  of 
launch  is  the  "final  time  for  the  launch  phase  and  the  Initial  time  for  the 
free-flight  phase. 

We  first  undertake  the  task  of  determining  the  orientation  of  the  Cu^u2U3 
frame  at  launch.  The  inertial  orientation  of  this  frame  is  defined  by  the 
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Euler  angles  <j>,  6 and  $ in  a 3-2-1  rotation  sequence.  Also,  Its  orientation 

relative  to  the  Pxyz  frame  is  as  a result  of  the  fact  that  it  is  the  centro idal, 

principal  frame  is  defined  by  the  two  small  angles  w2  and  ^3  which  must  be 

found  by  using  the  equations, 

L 

U2  - - y*o(zo+?-zc)(x-xc)dx/IT  (127a) 

0 

and 

L 

v*3  “ y*o(yo+n-yc)(x-xc)dx/IT  , (127b) 

0 

where  (x  , y ,z  ) is  the  location  of  the  rocket  center  of  mass  in  the  Pxyz 

O L < C 

reference  frame.  The  inertial  orientation  of  the  Pxyz  reference  frame  at  the 
instant  of  launch  is  defined  by  the  angles  ®L  and  $.  Thus,  the  inertial 
orientation  of  Cu^u^u^  reference  frame  at  the  instant  of  launch  can  be 
specified  using  either  0^,  $,  y2  and  y^,or  the  corresponding  values  of  <|>,  0 
and  <f>.  By  forming  the  rotation  matrices  which  are  associated  with  each  of 
these  sets  of  angles  and  equating  certain  elements,  the  following  equations 
for  ip,  0 and  d>  may  be  obtained: 

■ tan^Ku^0*  + U2s#)/[c0j+(iij  s*  “h2c<t)s0L]}  , (128a) 

0 * sin"1  [sQj+^d'-y^slOcOk]  (128b) 

and 


$ ■ tan_1{(-y,  s0  + s$c0  )/(-y  s©  + cQcO  )}.  (128c) 
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The  angular  velocity  of  the  Cu^u^  of  u are  related  to  0y  and  ^ 
through  the  approximate  equations. 


w.  - fl  , 

1 X 


W2  - y2  - P3Qx 


(129a) 

(129b) 


and 


a>3  - + ^2^  * 


which  are  approximate  only  in  that  m2  and  are  assumed  to  be  very  small 
angles  and  M2  and  M3  are  also  assumed  to  be  reasonably  small. 

Further,  the  location  of  the  center  of  mass  C with  respect  to  the  fixed 
OXYZ  frame  and  its  velocity  at  launch  are  needed.  Letting  Xc>  Yc  and  ZQ 
denote  the  inertial  coordinates  of  C, 


where 


+ x„  and  x*  in  the  x^-coordinate  of  point  P.  Moreover,  the 


XI‘C  " *Lp  + XC  and  \ 


Uj -components  of  the  velocity  of  the  center  of  mass  are,  up  to  launch. 


U‘\  * 


v " ~u3  u + yC  ‘ Vc 


(131a) 

(131b) 


and 


w - m2  u + + Vc* 


(131c) 


rx„i 

r cgt  0 s0. 1 

’1  0 0 ' 

"xT 

c 

L L 

Lc 

Yc 

- 

0 10 

0 c$  -s4> 

yc 

9 

-zc- 

.-S0L  0 ceL. 

_ 0 s$  c<t. 

-V 
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Finally,  Che  new  generalized  coordinates,  or  mode  shape  amplitudes,  and 
the  first  time  derivatives  can  be  found  by  using  the  following  equations: 


and 


where 


: " " u3b  2 + 22-ycbH  * 

(132a) 

! " y2b  2 + l 3 - zcb  l • 

(132b) 

: -^3b  5 + 1 i - ^Cb  i ■ 

(132c) 

! “ ^ ? + £ a - sa>  i * 

(132d) 

2b  ” E/XT  * 

(133a) 

3b  " £T  3/Xt  ’ 

(133b) 

o,  " ?T  2'Hj  ’ 

(133c) 

cb  “ ?T  3%  * 

(133d) 

L 

J(x-x  )£  dx  , 

(133e) 

o c 

T 

• * ^1  ^2  ' * * 

(133f ) 

L 

; - J $ dx 

(133g) 

and  the  elements  of  T are 
L 

t. . - j dx  (i,j-l,2,...,N), 

n 1 J 


- 3L-* 


(133h) 
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Equations  (132a)  and  (132b)  were  obtained  by  setting 

u2e  - - \^2e('x~xc)  + n (134a) 

u3e  " u2e^x“xC^  + 5 (134b) 

using  Eqs.  (104)  and  (124),  multiplying  each  of  Eqs.  (134)  by  4^  successively 

with  i - 1,2, . . . ,N  and  integrating  over  the  length  of  the  rocket.  In  Eqs. 

(134) , and  are  the  principal  axes  rotations  due  to  bending  only. 

Ze  3e 

Equations  (132c)  and  (132d)  were  obtained  from 

u2e  “ - u3e(x-xc)  + n (135a) 

and 

u3e  * P2e(x'xc^  + ? (135b) 

using  the  same  method  as  above. 

3.5  Comments  on  Simulation 

The  approximate,  ordinary  differential  equations  for  the  launch  phase 
and  the  free-fllght  phase  as  well  as  the  necessary  transformation  equations 
form  the  basis  for  a computer  program  which  was  developed  during  the  latter 
stages  of  this  study.  The  program  is  intended  to  be  used  in  studying  the 
effects  of  various  factors,  such  as  placement  of  the  supports  for  a rocket, 
on  the  motion  of  the  rocket.  It  was  hoped  that  this  program  would  be 
sufficiently  operational  in  time  to  compare  results  from  it  with  those  from 
the  two-body  rocket  model  program.  However,  the  program  is  not  yet  to  the 
stage  that  we  wish  to  publish  results  obtained  using  it. 
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From  the  preliminary  results  thus  far  obtained,  it  appears  that  the 
continuous  model  program  produces  numerical  values  of  the  right  orders  of 
magnitude  for  the  constrained  (on- launcher)  motions  of  a uniform  rocket 
model.  Also,  realistic  "launch"  values,  which  indicate  that  due  to  its 
bent  shape  at  launch  the  rocket's  principal  axes  have  nonzero  angular 
rates,  have  been  computed. 


SECTION  4.  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  Conclusions 


The  two-body  model  of  a flexible  free-fllght  rocket  yields  results  for 


angular  motion  of  the  nose  of  such  a rocket  which  agree  very  well  with  results 


derived  from  optical  lever  data.  To  achieve  agreement  in  the  case  considered, 


mechanical  thrust  misalignment  was  introduced  into  the  model  to  produce  oscil- 


lations in  the  time  histories  of  the  yaw  and  pitch  of  the  rocket  s nose  with  a 


frequency  close  to  the  spin  frequency.  Dynamic  mass  imbalance  of  an  equivalent 


amount  would  have  produced  the  same  results. 


The  very  good  agreement  between  theory  and  experiment  indicates  that  the 


model  is  both  qualitatively  and  quantitatively  valid.  This  in  turn  supports 


the  conclusion  that  the  transverse  bending  of  a free-flight  rocket  can  result 


in  significant  transverse  rates  at  launch.  As  an  example,  for  the  GEM #7  model, 


a one-milliradian  deflection  angle  of  the  nose  of  the  rocket  model  results  in 


a transverse  "steady-state"  rate  of  0.00816  rad/sec.  Furthermore,  a one- 


milliradian  per  second  angular  rate  of  the  nose  of  that  model  results  in  a 


transverse  "steady-state"  rate  of  0.000322  rad/sec.  When  consideration  is  given 


to  the  fact  that  the  angular  rates  due  to  bending  may  be  as  large  as  0.4  rad/sec. 


it  is  seen  that  "steady-state"  transverse  rates  of  free-flight  rockets  in 


excess  of  0.1  rad/sec  (=  100  mils/sec)  may  result  from  transverse  bending. 


The  magnitude  of  the  "steady-state"  transverse  rate  predicted  by  the  two- 


body  model  depends  significantly  upon  the  location  of  a rocket's  support  points 


(see  Table  2) , 


Although  the  presence  of  rather  severe  transverse  bending  of  free-flight 


rockets  has  been  documented12  the  results  obtained  from  the  two-body  model  for 
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on-launcher  morion  indicate  that  gravity  and  thrust  effects  alone  do  not  cause 
very  large  deflections  of  the  nose  of  the  free-flight  rocket  modeled,  hence, 
the  cause,  or  causes,  of  deflections  on  the  order  of  one-half  to  one  milli- 
radian,  which  have  been  observed,  are  still  not  known. 

As  far  as  the  continuous  rocket  model  is  concerned,  the  results  obtained 
thus  far  are  incomplete.  However,  the  use  of  two  different  types  of  mode 
shapes  appears  to  promise  significant  improvement  over  the  use  of  a single  type, 
because  the  transition  from  constrained  to  unconstrained  motion  is  modeled  in 
a much  more  realistic  manner. 

4.2  Recommendations 

The  following  recommendations  are  based  on  the  results  of  this  study. 

1.  Because  transverse  bending  of  a free-flight  rocket  can  result  in 
significant  transverse  rates  of  the  rocket  as  a whole  ("rigid  body"  rates) 
subsequent  to  launch,  a concerted  effort  should  be  made  to  determine  the  cause, 
or  causes,  of  bending  other  than  gravity  and  thrust  (as  it  was  modeled  herein). 

A by-product  of  such  an  effort  should  be  methods  for  reducing  transverse 
bending  and/or  designing  rockets  and  launchers  which  minimize  the  "rigid-body" 
rotation  rates  caused  thereby. 

I 1 

2.  The  two-body  model  computer  program  should  be  used  in  other  attempts 
to  duplicate  the  optical-lever-derived  pitch  and  yaw  time  histories  of  free- 
flight  rockets.  This  would  provide  additional  input  as  to  the  exactness  of  the 
model  and  also  lead  to  a better  understanding  of  the  problem  in  general. 

3.  The  continuous  rocket  model  program  should  be  used  in  the  same  manner 


as  the  two-body  one  with  the  primary  objective  of  comparing  the  results  obtained 
from  the  two  programs. 
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4.  Aerodynamic  forces  and  moments,  and  possibly  time  varying  inertia 
properties,  should  be  introduced  into  the  two-body  and  continuous  models 
used  for  the  free-f light  phase.  This  would  allow  simulation  of  the  rocket 
i from  launch  to  impact.  Hence,  the  way  in  which  transverse  bending  affects 

entire  trajectories  could  be  determined. 

! 


|j 
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APPENDIX  A 


DEFINITION  OF  THE  INTERNAL  FORCE 

The  internal  force  component  f(x)  in  the  rocket  is  assumed  to  be  due  to 
the  action  of  a force  with  the  magnitude  of  the  thrust  acting  at  the  head  end  of 
the  rocket  motor  as  shown  in  Fig.  A.l.  Since  no  aerodynamic  forces  are  considered  , 
the  only  other  forces  acting  on  the  rocket  are  the  force  due  to  gravity  and  a 
transverse  force  at  P which  is  due  to  the  turning  of  the  fluid  flow  from  the 
rocket  motor  by  a mechanically  misaligned  nozzle  and  a bent  rocket.  Letting 
denote  the  total  mass  of  the  rocket  and  0^  the  launcher  elevation  angle. 


•MTg  sin6l 


Fig.  A.l.  Simplified  model  for  determining  f(x). 


the  x-component  of  the  force  on  the  rocket  is  -Kpg  sin  0^  during  the  launch 
phase  and,  during  the  part  of  the  free-f light  phase  of  interest  here,  it  has 
essentially  this  value  because  the  pitch  and  yaw  angles  remain  small,  if  we 
do  not  consider  catastrophic  launches.  Also,  due  to  the  small  angular  rates 
and  transverse  velocity  components  the  acceleration  of  each  part  of  the  rocket 
in  the  x-direction  is  approximately  Ft/Hj,  ~ g sin  0^. 
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The  internal  force  at  a value  of  x less  than  x^  must  be  a tension  force 
and  so  is  positive  using  our  chosen  sign  convention.  This  force  must  be  of 
sufficient  magnitude  to  impart,  with  the  help  of  gravity,  to  the  part  of  the 
rocket  aft  of  its  location  the  acceleration  stated  above.  Thus,  since  the  mass 
accelerated  is 


Jo  dx. 


0<x<x  , 

- m 


+ f(x)  - m(x)  g sin  ©L  - m(x)  (F^,/^  ' 8 sin  ®L)  , 


f(x)  * (Ft/Mt)  m(x) . 


(A-l) 


(A-2) 


(A-3) 


At  x*x„  , the  force  f(x)  changes  to  one  of  compression.  At  x*x  it  must 
■f  “ M+ 

have  the  magnitude  -(F^/M^m^,  where  m^  is  the  mass  of  the  rocket's  nose 

section.  From  x*xM  to  x*L  the  magnitude  of  f(x)  decreases  because,  with 

increasing  x,  less  mass  is  to  be  accelerated.  Since 


“N  " / ° dx 
*M 

and  since  the  mass  accelerated  by  f(x)  for  x^<x<L  is 

x 

m(x)  * m^  - J a dx  , 

*M 

it  follows  that  for  x>x^  , 

x 

f (x)  - “(Ft/Mt)  [ffljj  *•  y o dx]  . 


(A-4) 


(A-5) 


(A-6) 
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For  the  special  case  of  constant  masses  per  length  for  the  motor  and 
for  the  nose. 


f(x)  - ~ CFt/Mt)  om  x , Oix<xM 
f(x)  - (Ft/H].)(L-x)on  , x^xxL. 


It  is  clear  from  Eqs.  (A-3)  and  (A-6)  that  ) 

- -F^/Hj,, while  f(0)  and  f(L)  are  zero.  Thus, 


/ 


3f/3x  dx  + 


3f/3x  dx 


(A-7a) 


(A- 7b) 

FT  VS  “d 


(A-8) 


The  force  f (x)  is  defined  as 

f(x)-  f(x)i+FT(  (03p  + az)j  +(e2p  - cty)k)6(x-0) 

where  6(o)  is  the  Dirac  delta  function  and  is  zero  for  all  x except  x-0  and 
is  1 at  that  value  of  x,  a and  a are  mechanical  thrust  misalignment  angles 

y z 

and  0_  and  0,  are  angles  of  rotation  of  the  aft  end  of  the  rocket  due  to 
2p  3p 

bending.  For  the  free-f light  phase,  the  unit  vector  i,  j and  fe  are  replaced 
by  u^»  u2  and  u^,  respectively. 
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MODE  SHAPES  FOR  THE  LAUNCH  PHASE 

The  mode  shapes  used  to  approximate  the  transverse  deflections  of  the 
continuous  rocket  model  during  the  launch  phase  are  those  for  a slender, 
uniform  beam  which  is  pinned  at  one  end,  point  P,  and  at  an  arbitrary  point 
Q as  shown  in  Fig.  B.l.  The  forward  end  is  free. 


Letting  El  denote  the  bending  stiffness  of  the  beam,  o its  mass  per  unit 
length,  and  L its  length,  the  equation  governing  its  transverse  free  vibration 


is 

ay*-  EIyIV  , (B-l) 

subject  to  the  boundary  conditions, 

y(0,t)  - 0 , (B-2a) 

y"(0,t)  - 0 , (B-2b) 

y(Xp,t)  ■ 0 , (B-2c) 

y"(L,t)  - 0 (B-2d) 

and 

y,M(L,t)*0  (B-2e) 

and  the  continuity  condition, 

■ y'Otq+.t)  • (®-3) 
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By  using  the  usual  method  of  separation  of  variables  whereby  we  set 


y - T(t)  Y(x) 


we  get  the  two  ordinary  differential  equations, 


(B-4) 


and 


T - -B**(EI/a)T 


YIV  - B4Y  , 


(B-5a) 


(B-5b) 


where  8 is  a separation  constant  which  must  of  course  have  only  certain 
values,  as  determined  by  the  boundary  conditions  (B-2)  and  the  continuity 
condition  (B-3) . 

The  general  solution  to  Eq.  (B-5b)  has  the  form. 


Y ■ A cosh  Bx  + B sinh  Bx  + C cos  Bx  + D sin  Bx 


(B-6) 


and  the  conditions  (B-2),  except  the  slope  continuity  conditions,  can  be 
used  to  obtain  two  functions  Y^  and  which  satisfy  those  conditions. 
Explicitly, 

Y^  ” sinx^l (sinBx) /sinBXq  ~(sinh8x)/sinh6xQ]  (B-7a) 

and 

Y2  - f(xQ)/h(XQ){[sin  6(1-5)  + sinh  8(1-5) ]/ [sin  Bo  + sinh  Bo] 

-[cos  B(l-x)  + cosh  8(l-x)]/[cos  So  + cosh  Bo]  , (B-7b) 


APPENDIX  C 


MODE  SHAPES  FOR  THE  FREE-FLIGHT  PHASE 


The  equations  which  govern  the  transverse,  vibrational  motion  of  a 
free-free,  slender  uniform  beam  is,  of  course,  the  same  as  that  given  in 
Appendix  B,  namely. 


ay  ■ Ely 


IV 


(C-l) 


where  again  o is  the  mass  per  unit  length,  y-y(x,t)  is  the  transverse 

displacement  at  x,  El  is  the  bending  stiffness,  (')*3(  )/3t  and 
IV  u u 

( ) “ 34(  )/3x  . However,  the  boundary  conditions  Eq.  (C-l)  are,  for 

a beam  of  length  L, 


y"’(0,t)  - 0 , 


yIV(o,t)  - o , 


y,M(L,t)  - o 


and 


IV 

y (L,t)  - 0 . 


(C-2a) 

(C-2b) 

(C-2c) 

(C-2d) 


By  setting  y - T(t)  Y(x),  and  separating  variables  in  Eq.  (C-l),  an 
equation  for  Y(x)  can  be  found.  Then,  by  using  the  boundary  conditions 
(C-2)  the  solution, 

Y(x)  ” (cos  6 - cosh  0) (sin  0x  + sinh  0x) 


-(sin  0 - sinh  0) (cos  0x  + cosh  0x) , 


(C-3) 
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1 


Ba  + sinh  Ba) (cos  8a  + cosh  Ba) 


The  continuity  of  the  slope  of  the  beam  at  x-Xp  requires  that  the 
eigenvalue  B satisfy  the  transcendental  equation. 


2(l+cos  Ba  cosh  Ba)/(sin  Ba  cosh  Ba  -cos  8a  sinh  Ba) 


denote  values  of  8 which  satisfy  Eq.  (B-9)  for  a 


given  value  of  a.  Then,  the  normalized  mode  shapes  i|>,  are  defined  by 


(B-10) 


(B-ll) 
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where  x - x/L,  8 « BL  and  Che  eigenvalue  g must  satisfy  Che  transcendental 
equation, 

cos  6 cosh  6 ” 1 , (C-4) 

may  be  obtained. 

We  require  here  only  the  bending  mode  shapes,  i.e.,  the  rigid  body  mode 
shapes  corresponding  to  g - 8^  * 0,  j*0,l,  are  not  used.  The  normalized  mode 
shapes  used  in  the  continuous  model  of  the  rocket  are 


♦ - Y(x,g  )/E  (j-1,2,...), 

where 


L 


0 


(C-5) 


(C-6) 
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TO:  Recipients  of  Subject  Report 

The  following  changes  should  be  made  to  subject  report: 

Page  4,  second  paragraph,  third  line:  "P"  should  be  "C^." 

Page  20,  Eq.  (38):  "I"  should  be  "I,." 
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Page  24,  first  line:  Delete  "negatives  of  the." 

Page  24,  Eq.  (48)  should  read:  x,  (t)  ■ E D E ^ x 

— h ■ - = — o 

Page  26,  Eq.  (55):  ”P(t)"  should  be  "P  (T)"  and  "2(t)" 

should  be  "2^(t)." 

Page  26,  Eq.  (56):  "P(t)"  should  be  "P^(t)m  and  and  the  [ J's 

on  the  (2,2)  element  of  P^(t)  should  be  deleted. 

Page  27,  second  line  from  bottom:  "sin  +"  should  be  "sin  4>  +." 

Page  28,  tenth  line:  "positive"  should  be  "negative." 

Page  28,  eighteenth  line:  "...apply  to..."  should  be  "...apply,  and..." 

and  ...quasi  rigid,..."  should  be  ". . .quasi-rigid , . . ." 

Page  28,  footnote:  Add  "and  their  displacements  relative  to  the 

main  part  of  the  spacecraft  are  small." 

Page  32,  seventh  line:  Second  "["  should  be  "(." 

ninth  and  tenth  lines:  Delete  "("  before  "d." 

thirteenth  line:  "P  " should  be  "0  ." 

o o 

fifteenth  line:  "0  " should  be  "T  ." 

o o 


Page  41,  Table  4:  Units  of  k should  be  "(Nt/m)." 
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dated  September  1978 

Page  71,  thirteenth  line:  Delete  first  "P"  within  brackets. 

P«rc  76,  third  line:  Delete  "is." 

Page  77,  first  line:  Change  to  read:  "The  Uj,  u^  and  u^  components 

of  w are  related  to  0^,  » *'3*  ^ >*3  as." 

sixth  line:  Add  equation  number  (129c). 

second  paragraph:  Add  equation  number  (130). 

second  paragraph,  seventh  line:  "in"  should  be  "is." 

Page  79:  Equations  (134a)  and  (134b)  should  read: 

U2c  - " - yCb 


“3,  ' - “2b(,-xC>  + ' " xCb  • 


Page  79,  sixth  line:  "p«  and  p,  " should  read  "p_.  and  p,.." 

2e  3c  2b  3b 

Page  79:  Equations  (135a)  and  (135b)  should  read: 

*2.  ' - S3b(,-xC)  + * - *Cfc 
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